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Nomenclature 

Ps 

time-averaged static pressure 

Cp 

dimensionless pressure coefficient 

Ps 

unsteady static pressure 

Cp s 

static pressure sensitive dimensionless 

Pt 

total pressure 


coefficient 

Pt 

time-averaged total pressure 

CpT 

target pressure sensitive dimensionless 

Pt 

unsteady total pressure 


coefficient 

Cpt 

total pressure sensitive dimensionless coefficient 

q 

dynamic pressure 

Cpta 

Re 

Reynolds number 

pressure coefficient 

Cptb 

U 

Downstream velocity 

pressure coefficient 


Cptc 

U 

time-averaged downstream velocity 

pressure coefficient 

Cpa 

U' 

unsteady downstream velocity 

pitch sensitive dimensionless pressure 


coefficient 

X 

wake survey downstream distance 

c P p 

yaw sensitive dimensionless pressure coefficient 

a 

flow angularity “pitch” 

d 

diameter of circular cylinder (wake survey) 

P 

flow angularity “yaw” 

i 

index variable (port number) 

y 

empirical coefficient (potential flow) 

K 

empirically determined coefficient 

e 

empirical coefficient (potential flow) 

P 

instantaneous mean pressure of off-axis pressure 

0 

included angle of incidence (potential flow) 


ports 

X 

probe port “cone angle” 

Pi 

pressure at the ith port of the seven-hole probe 

n 

dimensionless flow angularity coefficient 

Ps 

static pressure 

<t> 

probe port “clock angle” 


V 




Measurement of Air Flow Characteristics Using Seven-Hole Cone Probes 


TIMOTHY T. TAKAHASHI* 
Ames Research Center 


Summary 

The motivation for this work has been the development of 
a wake survey system. A seven-hole probe can measure 
the distribution of static pressure, total pressure, and flow 
angularity in a wind tunnel environment. The author 
describes the development of a simple, very efficient 
algorithm to compute flow properties from probe tip 
pressures. Its accuracy and applicability to unsteady, 
turbulent flow are discussed. 

Introduction 

For many years, multi-hole pressure probes have been 
used to measure flow angularity (refs. 1 and 2). The basic 
principle of operation is that a flow of given angularity 
and velocity at a given static pressure determines a unique 
pressure pattern across the various ports. The seven-hole 
probe has a specific construction advantage: six tubes of 
equal diameter fit exactly around a central tube of the 
same diameter. The geometry of the seven-hole probe 
is shown in figure 1. The resulting system is over- 
determined; there are more observations than states. 

The seven measured pressures determine the four output 
parameters: static pressure P s , total pressure Pt, and flow 
angularity a and p. 

In general, multi-hole pitot probes are used under the 
assumption that the flow of interest is steady. Great care 
is taken in the calibration process to ensure that the cali- 
bration flow is uniform and free from turbulence. How- 
ever, real flows, particularly wake flows, are far from 
steady. The response of the probe to unsteady flow must 
be documented. 

The decomposition algorithm, which reduces the seven 
pressures into the four air flow parameters, may be built 
using an ad hoc approach, an analytical flow model, or a 
mixture of theoretical and empirical rationale. This paper 
will describe the search for a procedure which combines 
good data quality with extreme processing efficiency. 


* This work was performed while the author held a National 
Research Counci l-N AS A Ames Research Associateship. 


Thanks to Pat Moriarty, Stanford University, for the hot- 
film anemometry data and to Tony Whitmore, Dryden 
Flight Research Center, for inspiration in the “triples” 
derivation. 

Application of Quasi-Steady Theory to 
Unsteady Flow 

There is a need to understand the applicability of any 
steady-state analysis, calibration, or computational 
procedure to the unsteady, highly turbulent flow often 
found in wake flow. In the past, many authors addressed 
this topic. A brief survey of some relevant publications is 
given here. Together, they describe issues important in 
the design of seven-hole probes. 

Goldstein (ref. 3) postulated that the response of a total 
head probe in a turbulent stream, Pt, is a measurement of 
mean total pressure P t plus a steady state contribution 
due to the unsteady velocity field: 

P t =P^ + P,' = P s l/2pU 2 + l/2pU' 2 

In other words, Goldstein considered static pressure a 
mean flow property. 

Jenkins (ref. 4) examined a free jet with a pitot tube 
equipped with a wide-bandwidth transducer. He corre- 
lated the unsteady pitot pressures with cross-wire hot-film 
velocimetry measurements. He found that the unsteady 
velocities as expressed by: 

U' = (U/2K)(P t '-P')/(Pt-^) 

were entirely representative of the hot- film signal. K is an 
empirically determined parameter; Jenkins reports typical 
values of 1.1 14 to 1.127. In other words, Jenkins consid- 
ered static pressure an instantaneous flow property. 

The relationship between the shape of a pitot probe tip 
and its time-averaged response in a turbulent flow was 
studied by Becker and Brown (ref. 5). Different probe 
geometries (tapered, round and square nosed) were found 
to have markedly different response to turbulence. This 
behavior may be used to advantage; an estimate of root- 
mean-square turbulence may be made by comparing the 



steady state stagnation pressures recorded by probes of 
differing geometry. 

Walshe and Gamer (ref. 6) studied the behavior of 
different probe configurations such as Pitot, Kiel, and 
five-hole probes in turbulent flow. They concluded that 
five-hole probes respond similarly to cowled pitot probes 
when measuring highly turbulent flow. Under these 
circumstances the mean dynamic pressure was judged to 
have up to 10% error. Flow angularity correlations were 
made comparing the five-hole probe used in a nulling and 
non-nulling configuration with the angularity of peak 
pressures obtained with an axial probe. Discrepancies 
were found between the three procedures when applied to 
turbulent flow. 

Whitmore (refs. 7 and 8) has extended the concept of 
multi-hole pressure probes to the construction of flush- 
mounted airdata systems. The HI-FADS system utilizes 
25 flush pressure ports mounted on an aircraft nosecone. 

It was calibrated for flight at both high angle of attack and 
large sideslip angles. At a 25 Hz computation rate, output 
from the HI-FADS system correlated well with flight 
data. One must note, however, that the flow over the 
aircraft nosecone is laminar. 

In summary, the existing literature implies the following: 
(1) when measuring a flow with frequency content 
beneath the pneumatic attenuation limit of the probe, 
quasi-steady techniques can measure both mean and 
unsteady flow properties; (2) above this frequency limit, 
the probes tend to overestimate the mean dynamic 
pressure; and (3) uncertainty exists when mean flow 
angularity measurements are made in turbulent flow. 

Despite these limitations, a probe calibrated under steady 
flow conditions remains a useful tool for measurement. 
Given sufficient frequency response, it can provide a 
metric of flow unsteadiness. A method of probe compu- 
tations consistent with steady state theory will be 
developed below. 

Potential Flow Model 

Whitmore (ref. 8) has demonstrated the applicability of 
a simple, hemispherical potential flow model to derive 
expressions for determining flow angularity. The pressure 
coefficient at the surface of a hemisphere is: 

C p (©) = 5/4 + 9/4cos 2 (0) (1) 


Cp(0) = e+ ycos 2 (0) (2) 

where £ and y are empirically determined for a given 
probe shape. 

The port pressures are: 

Pj =P s +q|cos 2 (©j) + esin 2 (0j)j (3) 


where 

0j = cos a cos P cos ^ +sinpsin^j sinfy 
+ sin a cos p sin X j cos 4> j 


(4) 


and 'ki and <|>i are coordinates of the pressure ports on the 
probe tip. 

Figures 2(a) and 2(b) show the response of three merid- 
ional probe tip pressure ports (ports 1, 7, and 4) as a 
function of pitch angle, a. In figure 2(a), the pressures are 
computed using equation 3 where the empirical coeffi- 
cient, £, is chosen so that the computed pressures at a = 0 
are consistent with experiment. Figure 2(b) shows the 
actual probe pressures measured using a 45° cone probe. 
A comparison reveals that the empirical potential flow 
model provides qualitative, but not quantitative, 
prediction of the port pressures. 


Method of Triples 

The potential flow model may be used to further examine 
the possibility of contriving dimensionless pressure ratios 
which are solely functions of flow angularity. The sim- 
plest pressure ratios, or “triples,” involve the pressures at 
three distinct ports: 

n ijk =( p i- p j) / (Pj-Pk) i*j (5) 

If one assumes that the pressures are governed by the 
potential flow model, equation 3, the pressures are: 

n ijk =(cos 2 ©j - cos 2 0j)/(cos 2 ©j - cos 2 © k j (6) 

Figures 3(a) and 3(b) demonstrate the differences between 
the actual probe response and simplified model given a 
pressure triple based upon three meridional pressure ports 
(ports 1, 7, and 4). Both the experimental data and the 
analytical model exhibit qualitative similarities including 
a singularity at a flow angularity of approximately one- 
half the probe cone angle. 


where 0 is the total flow incidence angle at the surface. 
Following Whitmore’s derivation, the pressures on quasi- 
hemispherical shapes, such as the cone probe or an 
aircraft nose cone, may be approximated by: 


The non-meridional triples are significantly less well 
behaved. As can be seen in figures 4(a) and 4(b), some of 
the non-meridional triples are multiply valued functions. 
n 12 3 = 0.25 may correspond to a flow angularity of -30°, 
-20°, -7°, and, possibly, +2°. 
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To reduce the system from an overdetermined state to one 
of one-to-one mapping, piecewise continuous pressure 
ratio functions must be used. Gallington (ref. 9) has found 
a complex, dimensionless pressure ratio is well behaved 
over a wide region of calibration space. A pair of pressure 
ratios, sensitive to pitch and yaw, are defined as functions 

of all seven pressures: 


c Pa7 = c Pta + ( c Ptb -C Ptc ) /2 

(7a) 

c Pp7 =(l/ V3)(Cp, b +Cpt c ) 

(7b) 

where 


Cp, a =(P4-Pl)/(P7-P) 

( 8 a) 

Cp,b=(P3-P6V(P7-P) 

( 8 b) 

Cptc=(P2-P5)/(P7-P) 

( 8 c) 

P = (l/ 6 XPi+P 2 +P 3 +P 4 + P5+P6) 

(9) 


The response of these coefficients to flow angularity is 
shown in figures 5 and 6 . These figures show that 
Gallington’s functions are well behaved and, more 
importantly, single valued over their entire calibration 
range. Provided that the flow over the entire probe tip is 
attached, all seven pressures are relevant. These coeffi- 
cients directly define the flow angularity. 

These calibration functions may be numerically inverted 
to define flow angularity, a and (}, in terms of the two 
coefficients. A pair of transfer functions, f a and fp, may 
be computed: 

a = f a( c Pai' C PPi) 

P = f p( c Pai* c PPi) 

They are tabulated in figures 7(a) and 7(b). 

The inversion process was accomplished using the 
following procedure. First, a subroutine was written to 
interpolate values of Cp a j and Cppj from the calibration 
data set given any arbitrary values of a and p. This 
subroutine was embedded in a minimization algorithm 
designed to find an a and p such that the interpolated 
values of Cp a j, Cppj match “target” values of these 
functions. In other words, to find the a and P corre- 
sponding to a given pair of values of the calibration 
functions: 


Step 1: Choose a target value for the flow angularity 

coefficients Cp a iT and Cppjx- 

Step 2 : Minimize: f(a,P) 

where 

f(a.P) = (Cp a j(a,P) - C PaiT ) 2 

2 

+ (Cppi (a, P) - CppiT ) 

Subject to: a,p bounded. 

Step 3: At local minimum? 

If YES, use derived a,p. 

If NO, let a,P go out of bounds 
(a = -999°, p = -999°). 

Step 4: Insert values of a,P into the 

fai(Cp<xi> Cppj) and fpi(C Pa j, Cppj) 
transfer function matrices. 

Step 5: Choose new values of Cp a iT and 

Cppj and i. 

Step 6: Go to Step 2. 

High Flow Angularity Coefficients 

For the high flow angularity situation, where the flow 
over the probe tip may have separated, unique pressure 
ratios must be formulated to avoid inclusion of ports on 
the lee side of the probe. Six pairs of coefficients may be 
formulated; each excludes specific neighboring ports. 
Because of the restricted domain of these functions, they 
need not be well defined about a = P = 0 . 

The rule set divides the calibration set into specific 
sectors. In Zilliac (ref. 10), the choice of sector is 
determined from an a priori inspection of the dominant 
pressure. If P 7 > Pj (i = 1 . . . 6 ), then Cp a 7 and Cpp 7 are 
utilized; if P 7 is not the dominant pressure, then the 
algorithm uses a separated flow pressure coefficient. Pi 
tends to exceed P 7 when the flow angularity reaches half 
of the cone angle (22.5° for a 45° cone probe). In reality, 
the flow will separate at much higher angles. This author 
concludes that Zilliac may have used separated flow 
coefficients prematurely. 
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A more lenient flow separation criterion may be devel- 
oped. Cp a 7 and Cpp7 are extremely well behaved over 
the entire calibration range. Consequently, even at high 
flow angularity Cp a and Cpp are a single valued function 
of flow angularity. The onset of flow separation may be 
identified by the angle found using Cp a 7 and Cpjft. 

For the case of a 45° cone probe, a cut-off point of 30° 
included angle was used. 

Gallington (ref. 9) also developed a set of secondary flow 
angularity coefficients which have the property of being 
bounded, ICp a i1 < 2 and ICppjl < 2, and free from singu- 
larities over their respective useful calibration ranges. 

Sector 1 (Pi dominant) 

Cp a i = n 126 7 = (Pi - P?)/[Pl - (P2 + P6)/2] (10a) 

Cppi = ni26=(P6 “ p 2V[Pl “ (P2 + P6V2] (10b) 

Sector 2 (P 2 dominant) 

Cpoc2 = ni 23 7 = (P2 - P?V[P2 - (Pi + P3)/2] (10c) 

Cpp 2 = n 123 = (Pi - P 3 )/[P 2 - (Pi + p 3 )/2] (lOd) 

Sector 3 (P 3 dominant) 

Cp a3 = n 234 7 = (P3 - P7V[ p 3 ” ( p 2 + p 4V2] (I0e) 

C P p 3 = n 234 = (P 2 - P4)/[P3 - (P2 + P4V2] 000 

Sector 4 (P4 dominant) 

Cpot4 = n 1237 = (P 4 - P 7 )/[P2 - (Pi + P3)/2] OOg) 
C P p4 = = (P3 - P5)/[P4 - (P3 + P5)/2] (10h) 

Sector 5 (P5 dominant) 

Cpo5 = IT4567 = (P5 - P7VIP5 - (P4 + P6)/2] (10i) 

Cpp5 = n 4 56 = (P4 - P6MP5 - (P4 + P6)/2] ( 10j) 

Sector 6 (P6 dominant) 

Cpa6 = ni567 = (P6 - P7>/[P6 - (?5 + Pl)/2] (10k) 

C P p 6 = nis6 = (P5 “ Pl)/[P6 ‘ (P5 + Pl)/2] 001) 

Figures 8(a) and 8(b) demonstrate the behavior of Cp a i 
for large positive a, and Cp a 4 for large negative a. Each 
of these functions is singular about a = 0 and shows 
strong asymptotic behavior for a > 30°. 


Algorithm— Part I: Flow Angularity 

The processing algorithm may be implemented as 
follows: 

Step 1: Compute the basic flow angularity 

coefficients: Cp a 7 and Cpp7. 

Step 2: If they are outside the range of the sector 7 

transfer function table, go to Step 5. 

Step 3: Perform bilinear interpolation on the 

sector 7 tables to determine a and (3. 

Step 4: Are a and p indicative of separated flow 

over the probe tip? 

If YES, continue on to Step S. 

If NO, then we are done. 

Step 5: Determine which pressure port has the 

largest positive pressure (the “dominant 
hole”)— this defines the sector, i. 

Step 6: Compute the appropriate coefficients: 

Cpai, Cppi. 

Step 7: Are these coefficients within the range of 

the tables? If not, then we have BAD data. 

Step 8: Perform bilinear interpolation on the 

separated flow tables to determine a and p. 

Computation of Static and Total Pressure 

With the flow angularity determined, the static and total 
pressure may be inferred from ratio of peak pressure to 
the mean of the pressures governed by unseparated flow. 
Static and total pressure coefficients are formulated in 
terms of the actual tunnel static pressure, P s , and total 


pressure, Pt: 

Cpsi = t(P2 + P6V2 - P s ]/[Pl - (P2 + P6V2] (Ha) 
Cpti = (Pi - P t y[Pl - (P 2 + P6)/2] (11b) 

Cp s 2 = [(Pi + P3V2 - P S ]/[P2 - (Pi + P3V2] (He) 
Cpt2 = (P2 - Pt)/[P2 - (Pi + P3)/2] (Hd) 

Cps3 = [(P 2 + P 4 )/2 - P S 1/[P3 - (P2 + P 4 y2] (He) 
Cpt3 = (P3 ~ P t y[P3 ~ (?2 + P 4 )/2] (HO 

Cp s4 = [(P3 + P 5 y2 - PsMP 4 - ( p 3 + P5)/2] (1 lg) 

Cpt 4 = (P 4 - P t y[P 4 - (P3 + P5)/2] (lib) 

Cps5 = t(P 4 + P6y2 - P&MP5 - (P4 + P6)/2] (Hi) 
Cpt5 = (P5 - Pt)/[P5 - (P 4 + P6)/2] ( 1 1 j) 
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Cps6 = [(P5+Pl)/2-P s ]/[P6-(P5 +Pl)/2J (Ilk) 

Cpt6 = (P6-PtV[P6-(P5+Pl)/2] (110 

and 

C Ps 7 =(P-P s y(P 7 -P) (llm) 

Cp t 7 =(P-P t y(P 7 -P) (1 In) 

In theory, the static and total pressure can be recon- 
structed for any flow angularity. For example, in 
unseparated flow the static and total pressures are 
reconstructed from the probe pressures and estimated 
flow angularity as: 

P t =P 7 -C Pt7 (a,p)(P 7 -P) (12) 

and 

P s =P 7 -Cp s7 (a,p)(P7 — P) (13) 

The reconstruction coefficients are well defined over 
their useful ranges; select matrices are tabulated in 
figures 9(a)-9(d). 

Special care must be taken when recording the data used 
for calibration to obtain a correct estimate of static 
pressure. While total pressure remains invariant in an 
inviscid flow, an actual wind tunnel will exhibit an axial 
static pressure gradient. If an incorrect estimate of static 
pressure is made at calibration time, the probes will 
estimate an incorrect dynamic pressure in operation. 

Algorithm — Part II: Static and Total 
Pressure 

In terms of the computational algorithm, recall that the 
flow angularity computations have determined both the 
appropriate sector, i, and the flow angularity (oc,p). These 
coefficients are used to perform a pair of bilinear interpo- 
lations upon the appropriate static and total pressure 
coefficient matrix. 

Step 9: Interpolate Cp s j from a,p and i. 

Step 10: Compute static pressure, P s . 

Step 11: Interpolate Cpxi from a,p and i. 

Step 12: Compute total pressure, P t . 

Implementation of this algorithm is extremely efficient. 
On a 486 PC, a computational throughput of several 
hundred reductions per second is realized. A real-time 
computation of flow properties is possible. 


Experimental Verification Turbulent 
Wake Flow 

Results from two recent experiments will be shown. Both 
are studies of the downstream wakes of circular cylinders. 
Figures 10(a) and 10(b) derive from an experiment made 
in a small developmental wind tunnel at the Ames Fluid 
Mechanics Laboratory; figures 1 1 and 12 derive from an 
experiment made at the Ames 7- by 10-Foot Subsonic 
Wind Tunnel No. 1. The natural pneumatic dissipation of 
the probe tubing limits the frequency response to approxi- 
mately 2 Hz. The Karman vortex street shed by the 
cylinders occurs at a much greater frequency: approxi- 
mately 250 Hz for the Re = 40,000 experiment and 
approximately 60 Hz for the Re = 400,000 experiment. 
The probe will consider the vortex street a source of 
unsteady flow. 

Figure 10(a) has the probe located 1.5 diameters down- 
stream of the cylinder. The two traces correspond to the 
total pressure in the free stream and the total pressure 
with the probe immersed to the side of the near-field 
vortex street. The probe is sampled first at 200 Hz, then at 
100 Hz; both in excess of the probe’s low-pass pneumatic 
limit. It can be seen that the probe gives a uniform 
response to either a steady state or turbulent flow. 

Figure 10(b) is of similar geometry, but with the probe 
moved farther back, to x/d = 5.5. The integration time is 
extended to 2.5 seconds, for an effective sample rate of 
0.4 Hz. The values for total pressure remain constant. 
These two figures imply that the probe exhibits essen- 
tially steady state behavior when immersed in an unsteady 
flow with a frequency content far above the probe limit. 

Additional data were taken to compare the performance 
of the seven-hole probe against a hot-film anemometer. 
For this experiment, the seven-hole probe was configured 
as it is to make two-dimensional wake surveys. The probe 
is slowly (0.5 inch per second) translated across the wake. 
Flow properties are measured in real time (three samples 
of 0.25 second integration time per second), and then 
gated to the desired spatial resolution. 

Figure 1 1 shows that there is close agreement between the 
seven-hole probe and hot-film anemometer in the free- 
stream velocity measurements. However, the seven-hole 
probe tends to consistently overestimate the velocity 
when immersed in the vortex street. There is some scatter 
in the seven-hole data, attributable to the short-time 
integration interacting with buffeting of the probe boom 
(unavoidable in production testing). Figure 12 reveals 
that the root-mean-square unsteady velocity levels in 
the vortex street are as high as 25% of the local mean 
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velocity. The overprediction of dynamic pressure, how- 
ever unwelcome, is consistent with Goldstein (ref. 3), 
Becker (ref. 5), and Walshe (ref. 6). 

Error Analysis 

A Monte Carlo simulation may be made to assess the 
sensitivity of the seven-hole probe algorithm to random 
error. A series of perturbations of increasing magnitude 
were statistically applied to basic pressure patterns 
corresponding to a flow angularity of a = 0°, (3 = 0° and 
a = 20°, p = -20°, respectively. The results are shown in 
table 1. A 1% deviation in port pressures tends to lead to 
an uncertainty in flow angularity of approximately 2° at 
low flow angles, but tends toward unpredictably large 
excursions at high flow angularity ; this is due to the 
shallow slopes of the separated-flow flow-angularity 
functions. Consequently, it is desirable to restrict the 
probe to as narrow a range of flow angularity as possible. 

Concluding Remarks 

Seven-hole cone probes are an accepted means to 
measure the essential mean properties of fluid flow. 

A rapid computational algorithm which incorporates as 
few as four bilinear interpolations and no more than six 
interpolations (16 to 24 array references and simple 
arithmetic) is presented. This has produced a computa- 
tional algorithm efficient enough to allow the flow 
properties to be implemented in real time. Unsteady flow 
properties with frequency content below the pneumatic 
limit of the probe may be resolved. Higher frequency 
unsteady flow will tend to bias the probe into indicating 
a dynamic pressure in excess of the actual value. Over 
the region of calibration space where the flow has not 
separated, the flow angularity computations are robust. 

At high flow angularity, the separated flow coefficients 
become very sensitive to small perturbations. The 
computed flow directionality at high flow angularity 
is best left for qualitative, rather than quantitative, 
presentation. 

Some issues require further investigation. There is a 
need to: 

• Develop design guidelines to define the optimal 
balance between probe size and frequency response. 

• Determine the limiting frequency where a quasi- 
steady calibration can be applied. 

• Find well behaved flow angularity coefficients which 
are more tolerant of random error, at high flow 
angularity, than the Gallington set. (This author has 
tried many with little success.) 


• Use seven-hole probes, in the absence of better flow 
angularity coefficients, in a semi-nulling configu- 
ration (where the peak flow angles are restricted). 

• Address the meaning of mean flow properties, 
particularly static pressure, in a turbulent flow. 

• Address the concept of unsteady, bandwidth limited 
static pressure, as reported by a real time multi-hole 
pitot probe, in turbulent flows. 
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Table 1. Monte Carlo error analysis — seven-hole probe 
computational algorithm 


Probe pressure 
port perturbation 

a = 0°/p = 0° 
RMS error 

a = 207p = -20° 
RMS error 

±0% 

±0° 

±0° 

±1/4% 

±0.5°/±0.5° 

±1.0°/±0.7° 

±1/2% 

±1°/±1° 

±2°/±1.5° 

±1% 

±2°/±2° 

±7.5°/±16° 



Figure 3(a). Simple, dimensionless pitch-sensitive 
pressure ratio as a function of pitch. Meridional “triple ” 
(n 174 = (Pi-P 7 )/(P 7 ~P 4 ))- Empirical (potential flow) model. 




Figure 3(b). Simple, dimensionless pitch-sensitive 
pressure ratio as a function of pitch. Meridional u triple " 
(n 174 = (Pi-P 7 )/(Pt-Pa))- Experimental results from 
45° cone probe. 


Figure 1. General schematic of a seven-hole probe. 



Figure 2(a). Pressure due to pitch. Empirical (potential 
flow) model. Meridional pressure ports (1, 7, and 4). 
e = 0; chosen to match experimental data at 0° pitch 
angle. 



Figure 4(a). Simple, dimensionless pitch-sensitive 
pressure ratio as a function of pitch. Non-meridional 
u triple ” (n 123 ). Experimental data from 45° cone probe. 
Note singularity near 0° pitch angle. 
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Figure 2(b). Pressure due to pitch. Experimental results. 
Meridional pressure ports (1, 7 , and 4)). Pressures from 
45° seven-hole cone probe calibration. 
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Figure 4(b). Simple, dimensionless pitch-sensitive 
pressure ratio as a function of pitch. Non-meridional 
“ triple w (n 672 ). Experimental data from 45° cone probe. 
Note singularities near -35° and +35° pitch angle. 
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Figured . Complex , dimensionless pitch-sensitive pressure 
ratio , Cp(x7, as a function of pitch (P = 0). Experimental 
results from 45° cone probe . 
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Figure 7(c). Pitch transfer function for separated flow. 
Sector 1. a = f a (Cp a i, Cppi). 



Figure 6. Complex, dimensionless yaw-sensitive pressure 
ratio, Cpp 7 , as a function of yaw (a = 0). Experimental 
results from 45° cone probe. 
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Figure 7(d). Yaw transfer function for separated flow. 
Sector 1. p = fp (Cpai, Cppi). 



Figure 7(a). Pitch transfer function for unseparated flow. 
Sector 7. a = f a (CpaJ, Cppy). 


Figure 8(a). Complex, dimensionless pitch-sensitive 
separated flow pressure ratio, Cpai , 3S 3 function of 
pitch (p = 0). 
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Figure 7(b). Yaw transfer function for unseparated flow. 
Sector 7. p = fp (Cp<x7> Cppj). 



Figure 8(b). Complex, dimensionless pitch-sensitive, 
separated flow pressure ratio, Cp<x 4 , as a function of 
pitch (p = 0). 
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Figure 9(a). Static pressure coefficients for separated flow. 
Sector 1. Cp s i. 
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Figure 9(b). Total pressure coefficients for separated flow. 
Sector 1. Cptf. 



Figure 10(a). Time history of flow. Data from cylinder 
wake test (Re = 40,000, x/d = 1.5). Two spatial locations 
chosen — one in the free stream, the other in the vortex 
wake. 200 and 100 Hz sampling rates. Note minor 
“ scatter * for 200 Hz sample rate (due to electrical noise) 
and minimal * scatter * for 100 Hz sample rate. There is no 
appreciable detection of the flow unsteadiness in the 
vortex wake. 
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Figure 9(c). Static pressure coefficients for unseparated 
flow. Sector 7. Cp s j. 
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Figure 9(d). Total pressure coefficients for unseparated 
flow. Sector 7. Cp\j. 
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Figure 10(b). Time history in wake. Data from cylinder 
wake test (Re = 40,000, x/d = 5.5). Two spatial locations 
chosen— one in the free stream, the other in the vortex 
wake. 2.5 second integration time, 0.4 Hz sample rate. 
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Figure 1 1. Axial velocity in cylinder wake (Re = 400,000, 
xJd = 7.5). Seven-hole and hot-film anemometry. 4 Hz 
data acquisition rate (0.1 second pressure integration 
time). 



Figure 12. Turbulence in cylinder wake (Re = 400,000, 
x/d = 7.5) from hot-film anemometry. 
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Appendix A — Additional Figures 
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Figure 13(a). Transfer function for separated flow. 
a = f a (Cp a l , Cppi). 
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Figure 13(b). Transfer function for separated flow. 

P = fp (Cpal < Cppi). 
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Figure 13(c). Transfer function for separated flow. 

a = fa {Cpo.2 ’ Cpp2)- 
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Figure 13(d). Transfer function for separated flow. 

P = fp(Cpa2 ■ Cpp2)- 
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Figure 13(e). Transfer function for separated flow, 
a -fa (Cpoc3 ■ Cpps). 
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Figure 13(f). Transfer function for separated flow. 

P = fp(Cpa2 . Cpp2). 



Cp<x4 


+ 

o 

o 

o 

1 

-17 -16 -16 -16 -16 

-16 

-16 

-16 

-16 

+ 0.25 

1 

-18 -18 -18 -18 -18 

-18 

-18 

-18 

-18 

+ 0.50 

1 

-21 -21 -21 -21 -21 

-21 

-21 

-20 

-20 

+ 0.75 

1 

-25 -24 -25 -25 -24 

-24 

-24 

-20 

-24 

+ 1.00 

1 

-30 -31 -31 -30 -30 

-30 

-30 

-30 

-30 

+ 1.25 


+45 -43 -42 -43 -43 

-41 

-43 

-41 

-40 

+ 1.50 

1 

+45 -45 -45 -45 -45 

-45 

-45 

-45 

+37 

+ 1.75 

1 

+45 -45 -45 -45 -45 

-45 

-45 

+45 

+39 

+ 2.00 

! 

+45 -45 -45 -45 -45 

-45 

-45 

+45 

+41 


- 2.0 - 1.00 + 0.0 + 1.0 + 2.0 

Cp|34 

Figure 13(g). Transfer function for separated flow. 

a = f a (Cpa4 . Cpp4). 
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Figure 13(h). Transfer function for separated flow. 

P = fp(Cpa4 . Cpp 4 ). 
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Figure 13(i). Transfer function for separated flow. 

a = fct (Cpa5 , Cpps). 
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Figure 13(j). Transfer function for separated flow. 

P = fp(Cpa5 ■ Cpps)- 
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Figure 13(k). Transfer function for separated flow. 
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Figure 131(1). Transfer function for separated flow. 

P = fp(Cpa6, Cpp6). 
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Figure 14(a). Static pressure coefficients, Cp S f. 
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Figure 14(b). Total pressure coefficients, Cpti. 
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Figure 14(c). Static pressure coefficients, sector 2. 
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Figure 14(d). Total pressure coefficients, sector 3. 
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Figure 14(e). Static pressure coefficients, sector 3. 
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Figure 14(f). Total pressure coefficients, sector 3. 
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Figure 14(g). Static pressure coefficients, sector 4. 
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Figure 14(h). Total pressure coefficients, sector 5. 
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Figure 14(i). Static pressure coefficients, sectors. 
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Figure 14(j). Total pressure coefficients, sector 5. 
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Figure 14(k). Static pressure coefficients, sector 6. 


P 


+ 5.0 

| - 0.1 

- 0.0 

- 0.1 

- 0.6 

- 4.9 

+ 15.0 

- 0.1 

+ 0.0 

- 0.0 

- 0.4 

- 3.4 

+ 25.0 

- 0.1 

+ 0.0 

- 0.0 

- 0.5 

- 4.6 

+ 35.0 

j - 0.1 

- 0.1 

- 0.2 

- 1.0 - 11.4 

+ 45.0 

| - 0.2 

- 0.2 

- 0.5 

- 1.9 + 86.8 


+ 

- 45.0 

- 35.0 

- 25.0 

- 15.0 

- 5.0 


a 


Figure 141(1). Total pressure coefficients, sector 6. 


14 



Appendix B — Calibration Code (Microsoft QuickBasic) 

DECLARE FUNCTION min! (A!, B!, C! ( d!, e!, F!, g!) 

REM *************************************************************** 
REM *************************************************************** 
REM ** 

REM ** SHPCAL.BAS - program to build a calibration file for SHP.BAS 
REM ** 

REM ** Program by Timothy Takahashi - National Research Council 
REM ** 

REM ** Last Revision : October 5th, 1995 
REM ** 

REM ** PC-Compatible / MS-DOS 6.2 / Microsoft Quick Basic v. 4.5 
REM ** 

REM *************************************************************** 
REM *************************************************************** 

REM $ DYNAMIC 

REM >>> Dimension working arrays 


OPTION BASE 1 


DIM 

tablel (7 , 

20, 

20) 

DIM 

table2 (7, 

20, 

20) 

DIM 

table3 { 7 , 

20, 

20) 

DIM 

table4 (7, 

20, 

20) 

DIM 

pres (400, 

8) 



DIM alpha (4 00), beta (400), NDpres(8) 

REM »> Open output files 

OPEN "O" , 2, 11 c : \qbasic\pcdas\shpcal ,dat " 

REM >>> Initialize working arrays 

FOR probe% = 1 TO 3 
PRINT " PROBE 11 ; probe% 

PRINT "Initializing. . . " 

FOR i = 1 TO 20 


FOR j = 1 

TO 

20 

FOR k = 1 

TO 

7 

tablel (k, 

3 , 

i) = 99 

table2 (k, 

j / 

i) = 99 

table3 (k, 

j , 

i) = 99 

table4 ( k, 
NEXT k 

j , 

i) = 99 


NEXT j 
NEXT i 


REM >>> Read calibration file into working arrays 


PRINT "READING CALIBRATION PRESSURES" 

OPEN "I", 1, "c:\qbasic\pcdas\calpr.dat" 

i% = 1 
110 : 

IF EOF ( 1 ) THEN 120 
INPUT #1, alpha (i%) 

INPUT #1, beta(i%) 

INPUT #1, pres ( i% , 8) 

FOR P% = 1 TO 3 
IF P% = probe% THEN 
INPUT #1, pres (i%, 1) 

INPUT #1, pres (i% , 2) 

INPUT #1, pres ( i% , 3) 

INPUT #1, pres ( i% , 4) 

INPUT #1, pres ( i% , 5) 

INPUT #1, pres ( i% , 6) 

INPUT #1, pres ( i% , 7) 

ELSE 

INPUT #1, x 
INPUT #1, x 
INPUT #1, x 
INPUT #1, x 
INPUT #1, x 
INPUT #1, x 
INPUT #1, x 
END IF 
NEXT P% 
i% = i% + 1 
GOTO 110 

120 : ncal = i% - 1 

CLOSE 1 

PRINT "Searching for Reference Total Pressure" 

REM >>> Determine the reference total and static pressures. 

REM »> The reference total is p7/p8 at alpha=beta=0 . 

FOR i = 1 TO ncal 

IF (ABS (alpha (i ) ) < 1 AND ABS(beta(i)) < D THEN 

IER = 1 

PTOT = pres ( i , 7) / pres(i, 8) 

PDEN = pres(i, 8) 

PBARL = (pres(i, 1) + pres(i, 2) + pres(i, 3) + pres(i, 4) + pres(i, 

5 ) + pres ( i , 6 ) ) / 6 ! 

GOTO 90 
END IF 

NEXT i 

PRINT "No reference total pressure found" 
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STOP 


90 : 

PRINT "REFERENCE TOTAL PRESSURE = PTOT * PDEN; "PSI" 

PST = 0 

PSTAT = PST / PDEN 

REM >>> Compute Calibration Coefficients 
PRINT "Computing Coefficients" 

FOR i = 1 TO ncal 

REM >>> non-dimensionalize pressures 
FOR j = 1 TO 8 

NDpres(j) = presfi, j) / pres(i, 8) 

NEXT j 

ix% = (alpha(i)) / 5 + 10 
iy% = (beta(i)) / 5 + 10 

REM >>> Check to make sure pressure at hole 7 is not the lowest 
REM >>> pressure (Second check for separation at hole 7). 

PMIN = min (NDpres (1 ) , NDpres (2), NDpres{3), NDpres(4) / NDpres (5), 
NDpres ( 6 ) , NDpres ( 7 ) ) 

IF (NDpres (7) = PMIN) THEN 150 

REM >>> Compute the Coefficients for each Sector 

REM >>> Zone 1. 

PBAR1 = (NDpres (2) + NDpres (6)) / 2! 

PDIFF = NDpres (1) - (NDpres (2) + NDpres (6)) / 2! 

IF ( ABS ( PDIFF ) <= .0000001#) THEN PDIFF = .0000001# 

CT1 = (NDpres (1) - NDpres (7)) / PDIFF 
CFl = (NDpres ( 6 ) - NDpres{2)) / PDIFF 
CPTOT1 = (NDpres (1) - PTOT) / PDIFF 
CPSTA1 = ( PBAR1 - PSTAT) / PDIFF 
tableKl, ix%, iy%) = CT1 
table2 { 1 , ix%, iy% ) = CFl 
table3(l, ix%, iy%) = CPTOT1 
table4 ( 1 , ix%, iy%) = CPSTA1 

REM >>> Zone 2 . 

PDIFF = NDpres (2) - (NDpres (1) + NDpres (3)) / 2! 

PBAR2 = (NDpres (1) + NDpres (3)) / 2! 

IF (ABS (PDIFF) <= .0000001#) THEN PDIFF = .0000001# 

CT2 = (NDpres (2) - NDpres (7)) / PDIFF 
CF2 = (NDpres (1) - NDpres (3)) / PDIFF 
CPTOT2 = (NDpres (2) - PTOT) / PDIFF 
CPSTA2 = (PBAR2 - PSTAT) / PDIFF 
tablel(2 / ix% , iy%) = CT2 
table2(2 / ix% , iy% ) = CF2 
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table3 (2 , ix%, iy%) = CPT0T2 
table4 (2 r ix%, iy%) = CPSTA2 


REM »> Zone 3 . 

PDIFF = NDpres {3 ) - (NDpres(2) + NDpres(4)) / 2! 
PBAR3 = (NDpres (2) + NDpres (4)) / 2! 

IF (ABS (PDIFF) <= .0000001#) THEN PDIFF = .0000001# 
CT3 = (NDpres (3) - NDpres (7)) / PDIFF 
CF3 = (NDpres (2) - NDpres (4)) / PDIFF 
CPTOT3 = (NDpres (3) - PTOT) / PDIFF 
CPSTA3 = (PBAR3 - PSTAT) / PDIFF 
tablel(3, ix%, iy%) = CT3 
table2(3, ix%, iy%) = CF3 
table3 ( 3 , ix%, iy%) = CPTOT3 
table4 ( 3 , ix%, iy%) = CPSTA3 


REM »> Zone 4 . 

PDIFF = NDpres (4) - (NDpres (3) + NDpres (5)) / 2! 
PBAR4 = (NDpres (3) + NDpres (5)) / 2! 

IF (ABS (PDIFF) <= .0000001#) THEN PDIFF = .0000001# 

CT4 = (NDpres (4) - NDpres (7)) / PDIFF 

CF4 = (NDpres (3) - NDpres (5)) / PDIFF 

C PTOT 4 = (NDpres (4) - PTOT) / PDIFF 

CPSTA4 = ( PBAR4 - PSTAT) / PDIFF 

tablel(4, ix%, iy%) = CT4 

table2(4, ix%, iy%) = CF4 

table3 (4, ix%, iy%) = C PTOT 4 

table4 ( 4 , ix%, iy%) = CPSTA4 


REM »> Zone 5 . 

PDIFF = NDpres (5) - (NDpres (4) + NDpres (6)) / 2! 
PBAR5 = (NDpres (4) + NDpres (6)) / 2! 

IF (ABS ( PDIFF) <= .0000001#) THEN PDIFF = .0000001# 

CT5 = (NDpres (5) - NDpres (7)) / PDIFF 

CF5 = (NDpres (4) - NDpres (6)) / PDIFF 

C PTOT 5 = (NDpres (5) - PTOT) / PDIFF 

CPSTA5 = ( PBAR5 - PSTAT) / PDIFF 

tablel(5, ix%, iy%) = CT5 

table2(5, ix%, iy%) = CF5 

table3 ( 5 , ix%, iy%) = CPTOT5 

table4 (5, ix%, iy%) = CPSTA5 


REM »> Zone 6 . 

PDIFF = NDpres (6) - (NDpres (5) + NDpres (1)) / 2! 
PBAR6 = (NDpres (5) + NDpres (1)) / 2! 

IF (ABS (PDIFF) <= .0000001#) THEN PDIFF = .0000001# 

CT6 = (NDpres (6) - NDpres (7)) / PDIFF 

CF6 = (NDpres (5) - NDpres (1)) / PDIFF 

CPTOT6 = (NDpres (6) - PTOT) / PDIFF 

CPSTA6 = (PBAR6 - PSTAT) / PDIFF 

tablel(6, ix%, iy%) = CT6 

table2(6, ix%, iy%) = CF6 

table3 ( 6 , ix%, iy%) = CPT0T6 

table4(6, ix%, iy%) = CPSTA6 
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REM >>> Zone 7 . 

PBARL = (NDpres(l) + NDpres(2) + NDpres(3) + NDpres(4) + NDpres(5) + 
NDpres (6) ) / 6 ! 

PDIFFL = NDpres (7) - PBARL 
CA1 = (NDpres (4) - NDpres (1)) / PDIFFL 

CA2 = (NDpres (3) - NDpres (6)) / PDIFFL 

CA3 = (NDpres (2) - NDpres (5)) / PDIFFL 

CBL = 1! / SQR ( 3 ! ) * (CA2 + CA3 ) 

CAL = CA1 + (CA2 - CA3 ) / 2! 

PDIFF = PDIFFL 
PBAR7 = PBARL 

IF ( ABS ( PDIFF ) <= .0000001#) THEN PDIFF = .0000001# 

CT7 = CAL 

CF7 = CBL 

cptot7 = (NDpres (7) - PTOT) / PDIFF 
cpsta7 = ( PBAR7 - PSTAT) / PDIFF 
tablel(7, ix%, iy%) = CT7 

table2(7, ix%, iy%) = CF7 

table3(7, ix%, iy%) = cptot7 

table4(7, ix% , iy%) = cpsta7 


150 NEXT i 

PRINT "EVALUATING STATIC AND TOTAL PRESSURE COEFFICIENTS" 

REM >>> Write Static and Total Pressure Coefficient Arrays to Disk 

FOR ix% = 1 TO 20 
FOR jx% = 1 TO 20 
FOR ndhol% = 1 TO 7 

WRITE #2, probe%, ndhol%, ix%, jx%, table3 (ndhol% , ix%, jx%), table4 (ndhol% , ix%, 

jx%) 

NEXT ndhol% 

NEXT jx% 

NEXT ix% 


REM >>> Invert flow angularity matrix 
FOR ndhol% = 7 TO 1 STEP -1 

PRINT "INVERTING CALIBRATION MATRIX FOR SECTOR ndhol% 
REM >>> Different ranges for different sectors 
IF ndhol% = 7 THEN 


xO = -4.5: xl = 

4.5: 

x2 = 

.5 

x3 = -4.5: x4 - 
ELSE 

4.5: 

x5 = 

.5 

xO = 0 : xl = 2 : 

x2 = 

.25 


x3 = -2: x4 = 2; 
END IF 

: x5 

= .5 



FOR CpAlpha = xO TO xl STEP x2 
FOR CpBeta = x3 TO x4 STEP x5 
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ix% = 0 
jx% = 0 
eO = 9999999 

REM >» Search for minimum error in three phases 

REM >» Phase 1 - simple search w/o interpolation 

FOR i% = 1 TO 19 
FOR j% = 1 TO 19 

e = (tablel (ndhol% , i%, j%) - CpAlpha) A 2 + (table2 (ndhol% , i%, j%) - 
CpBeta) A 2 

IF e < eO THEN ix% = i%: jx% = j%: eO = e 
NEXT j% 

NEXT i% 

IF (ix% = 0) OR (jx% = 0) THEN BEEP: alpha = -999: beta = -999: GOTO 1001 

IF ( ix% = 1) THEN ix% = 2 
IF (ix% = 19) THEN ix% = 18 
IF (jx% = 1) THEN jx% = 2 
IF (jx% = 19) THEN jx% = 18 

REM »> Phase 2 - interpolate 


eO = 99999 


FOR il = ix% - 1 TO ix% + 1 STEP .25 
FOR jl = jx% - 1 TO jx% + 1 STEP .25 
iy% = INT(il): iz = il - iy% 
jy% = INT(jl): jz = jl - jy% 


A = tablel (ndhol%, iy% + 1, jy% + 1) * iz 

- iz) * jz + tablel (ndhol% , iy%, jy%) * d 
1, jy%) * iz * (1 - jz) 

B = table2 (ndhol% , iy% + 1, jy% + 1) * iz 

- iz) * jz + table2 (ndhol% , iy%, jy%) * d 
1, jy%) * iz * (1 - jz) 


* jz + tablel (ndhol%, iy%, jy% +1) * d 

iz) * (1 - jz) + tablel (ndhol%, iy% + 

* jz + table2 (ndhol% , iy%, jy% + 1) * d 

iz) * (1 - jz) + table2 (ndhol%, iy% + 


e = (A - CpAlpha) ~ 2 + (B - CpBeta) ~ 2 
IF e < eO THEN i 2 = il: j2 = jl: eO = e 


NEXT jl 
NEXT il 

eO = 99999 

IF (i2 < 1.25) THEN i2 = 1.25 
IF ( i2 > 18.75) THEN i2 = 18.75 
IF { j 2 < 1.25) THEN j2 = 1.25 
IF ( j 2 > 18.75) THEN j2 = 18.75 

REM »> Phase 3 - Interpolate over a finer grid 

FOR i3 = i2 - .25 TO i2 + .25 STEP .0625 
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FOR j 3 = j 2 — .25 TO j2 + .25 STEP .0625 

iy% = INT ( i 3 ) : iz = i3 - iy% 

jy% = INT ( j 3 ) : jz = j3 - jy% 

A = tablel (ndhol% , iy% + 1, jy% +1) *iz*jz+ tablel (ndhol% , iy%, jy% +1) * (1 

- iz) * jz + tablel (ndhol%, iy%, jy%) * (1-iz) * (1-jz) + tablel (ndhol% , iy% + 

1, jy%) * iz * (1 - jz) 

B = table2 (ndhol% , iy% + 1, jy% +l)*iz*jz+ table2 (ndhol%, iy% , jy% +1) * (1 

- iz) * jz + table2 (ndhol%, iy%, jy%) * (1 - iz) * (1 - jz) + table2 (ndhol% , iy% + 

1- jy%) * iz * (1 - jz) 

e = (A - CpAlpha) A 2 + (B - CpBeta) A 2 

IF e < eO THEN i4 = i3 : j 4 = j 3 : eO = e 

NEXT j 3 
NEXT i3 

alpha = (i4 - 10) * 5 
beta = (j4 - 10) * 5 

PRINT USING "SECTOR # | CPALPHA +#.### CPBETA +#.### | ALPHA +##.# BETA +##.#"; 
ndhol%; CpAlpha; CpBeta; alpha; beta 

IF eO > 1 THEN BEEP: alpha = -999: beta = -999: GOTO 1001 

1001 : 

WRITE #2, probe%, ndhol%, CpAlpha, CpBeta, alpha, beta 
NEXT CpBeta 
NEXT CpAlpha 

NEXT ndhol% 

PRINT 

NEXT probe% 

CLOSE 

END 

REM $ STATIC 

FUNCTION min (A, B, C, d, e, F, g) 
n = 0 

IF A < B THEN n = A ELSE n = B 

IF C < n THEN n = C 

IF d < n THEN n = d 

IF e < n THEN n = e 

IF F < n THEN n = F 

IF g < n THEN n = g 

min = n 

END FUNCTION 
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Apendix C — Data Acquisition Code (Microsoft QuickBasic) 

DECLARE SUB readanalog { ) 

DECLARE SUB readscales ( ) 

DECLARE SUB analoginit ( ) 

DECLARE SUB ReadCal () 

DECLARE SUB getdata { ) 

DECLARE SUB DebugProbes { ) 

DECLARE SUB crunch <P!(), 0!{), probe%) 

DECLARE FUNCTION max! (a!, b! , c\ , d! , el, f ! , gi) 

DECLARE FUNCTION min! (a!, b!, c!, d!, e!, f|, g!) 

DECLARE SUB DecodeBuffer (buffer$) 

DECLARE SUB Talk8400 (CMD$) 

DECLARE SUB SloSurvey ( ) 

DECLARE SUB survey ( ) 

DECLARE SUB Main { ) 

DECLARE SUB PositionMenu ( ) 

DECLARE SUB SmartMove () 

DECLARE SUB StillMoving () 

DECLARE SUB Init () 

DECLARE SUB SetUpCompumotor ( ) 

DECLARE FUNCTION keystroke% (0$) 

DECLARE SUB Menu (m$) 

DECLARE SUB INIT1 ( ) 

DECLARE SUB TalkCompumotor (Hstep&, Vstep&, Lstep&) 

DECLARE SUB GetPosition ( ) 

DECLARE SUB move to (Hnew! , Vnew! , Lnew! ) 


' Common GPIB status variables 

COMMON SHARED /NISTATBLK/ ibsta%, iberr% , IBCNT%, IBCNTL& 

’ GPIB status bit vector : 

status variable ibsta and wait mask 

CONST EERR = &H8000 1 Error detected 

CONST TIMO = &H4000 ’ Timeout 

CONST RQS = &H800 ’ Device requesting service 


' GPIB Subroutine Declarations 

. DECLARE SUB IBDEV ( BYVAL BDID% , BYVAL PAD% , BYVAL SAD%, BYVAL TMO% , BYVAL EOS% 
BYVAL EOT% , ud% ) 

DECLARE SUB IBCLR (BYVAL BD%) 

DECLARE SUB IBRDF (BYVAL BD% , FLNAME$) 

DECLARE SUB IBRD (BYVAL BD% , RD$ ) 

DECLARE SUB ibrsp (BYVAL BD% , spr%) 

DECLARE SUB ibwait (BYVAL BD% , BYVAL MASK%) 

DECLARE SUB IBWRT (BYVAL BD% , WRT$ ) 
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' Common block and Subroutine Declarations required 
' by the National Instruments 488 interface card (PC2A) 


DIM SHARED psi%, spx% , echo% 

DIM SHARED Nvals, PktInfo(19), pktmsg$, pktXdat (4096) 
DIM SHARED buffer3$, x$ 

DIM SHARED volts(16), scale(lO) 


REM >>> Sector 7 angularity matrix 
DIM SHARED tablela(2, 19, 19) 

DIM SHARED table2a(2, 19, 19) 

REM »> Sector 1-6 angularity matrix 
DIM SHARED tablelb(2, 6, 9, 9) 

DIM SHARED table2b(2, 6, 9, 9) 

REM »> Total/Static Pressure matrix 
DIM SHARED table3(2, 6, 20, 20) 

DIM SHARED table4(2, 6, 20, 20) 


buf fer3$ = SPACE${4096) 


' Execuatable code begins here. 


DIM SHARED xpos 
xpos = -999 

DIM SHARED pt$ , test$ 

DIM SHARED tX, ty, Tz 

CONST delay = 149 

REM * SLOW- SURVEY. BAS - SURVEY RIG CONTROL MODULE w/ 7 -HP DATA ACQ. 

t 

' Serial port mnn 

CONST portl$ = "COM2 :2400,N, 8, 1, DS0 , RS , CD0 " 1 compumotor 3000 link 

CONST motor% =0 ' = 1 to disable motor control, = 0 to enable 

•Set-up’s for MetraByte Parallel Ports 
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• NOTE : &H90 = PA-Input , PB-Output , PC ( 0-3 ) -Output , PC ( 4-7 ) -Output 
' Metrabyte to encoder wiring diagram: 

' PA (0-7) input < D ( 0-7 ) output from each encoder 

• pcO output > AO (pin #4) of each encoder 

' PCI output > A1 (pin #3) of each encoder 

' PC 2 output > Strobe (pin #6) of each encoder 

■ PC (3-7) > Device Select (pin #5) for encoder 1-5 

CONST BASEADDR = &H2E0 ' MetraByte Base Address (from switch settings 

[ 00111111 ] ) 

CONST PA = BASEADDR + 0 
CONST PB = BASEADDR + 1 
CONST PC = BASEADDR + 2 

CONST CTRL = BASEADDR +3 'MB Port Addresses 

CONST CONTROL = &H90 'MetraByte Port-configuration Contol Word 

'General Rig Parameters 

CONST Hscale = 1000 'Scale factor to turn Inches into Compumotor 3000 counts 

CONST VScale = 1000 'Scale factor ... 

CONST LScale = 1 / 3 * 1000 'Scale factor ... 


CONST WELOC = 15 
CONST LVELOC =10 

CONST HAccel =10 
CONST VAccel =10 
CONST LAccel = 10 

' SURBEY RIG HOME POSITION 

CONST HHome = 0 'Horizontal Home Position in inches 

CONST VHome = 0 'Vertical Home Position in inches 

CONST LHome = 0 'Longitudinal Home Position in inches 

' SURVEY RIG MAXIMUM EXTENSION POSITION (inches) 

CONST HLoLim = -32 'Horizontal Low Software Limit in Inches 
CONST HHiLim = 32 'Horizontal High Software Limit in Inches 

CONST VLoLim = -1 'Vertical Low Software Limit in Inches 

CONST VHiLim = 60 'Vertical High Software Limit in Inches 

CONST LLoLim = -16 'Longitudinal Low Software Limit in Inches 

CONST LHiLim = 32 'Longitudinal High Software Limit in Inches 

DIM SHARED closedloop% 


CONST false = -1 
CONST true = 1 

DIM SHARED Vdir, ERP, hveloc 
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DIM SHARED P$(7), ACT (4), ESB$(4, 8) 


encode control arrays 


CALL Init 
CALL ReadCal 
CALL SetUpCompumotor 
CALL Main 


CHDIR M c : \qbasic\pcdas " 

RUN "MENU. BAS” 

SUB analoginit 

DIM param% ( 60 ) , DAT%(32000) / gain%(15), nn ( 15), pp(15), offset(15) 


param%{0) = 0 ' Board number 

param%(l) = &H300 1 Base I/O address 

param% ( 4 ) = 2 1 IRQ level : IRQ2 

param%(6) = 20: param%(5) = 20: xmax = 2000 ’ sample at 500hz 

param%(7) = 0 1 Trigger mode, 0 : pacer trigger 

param%(8) = 0 ' Non-cyclic 

param%(10) = VARPTR ( DAT% ( 0 ) ) ' Offset of A/D data buffer A 

param%(ll) = VARSEG ( DAT% ( 0 ) } ' Segment of A/D data buffer A 

param%(12) = 0 ' Data buffer B address, if not used, 

param%(13) = 0 ' must set to 0. 

param%(14) = 4096 
param%(15) = 0 
param%(16) = 15 

param%(17) = &HFF ' gain table 

param%(18) = VARPTR (gain% (0) ) ' Offset of gain table buffer 

param%(19) = VARSEG (gain% ( 0 ) } ' Segment of gain table buffer 


gain% ( 0 ) = 

0: 

pp ( 0 ) = 

10: offset(0) = 

-5: 

+/- 5v external input 

gain%(l) = 

0: 

pp(l) = 

10: offset(l) = 

-5: 

Pa 

gain%(2) = 

8: 

pp ( 2 ) = 

20: offset(2) = 

-10: 

' Pr-Pa 

gain%(3) = 

8: 

PP ( 3 ) = 

20: offset{3) = 

-10: 

' Pr-Ps 

gain%(4) = 

4: 

pp ( 4 ) = 

10: offset ( 4 ) = 

0: ' 

turntable 

gain% ( 5 ) = 

0: 

PP ( 5 ) = 

10: offset ( 5 ) = 

-5: 

inclinometer 

gain%(6) = 

6: 

pp ( 6 ) = 

.1: offset(6) = 

0: 1 

PTplenum barocel 

gain%(7) = 

6: 

PP (7 ) = 

.1: offset{7) = 

0: ' 

Tplenum thermocouple 

gain%(8) = 

4: 

pp(8) = 

10: offset ( 8 ) = 

0: ' 

CFM plenum 

gain%(9) = 

9: 

pp ( 9 ) = 

2: offset (9) = - 

-1: 1 

1 PS forebody 

gain% { 10 ) 

= 9 

pp(10) 

= 2: offset(10) 

= -1 


gain% (11) 

= 9 

PPdl) 

= 2: offset ( 11 ) 

= -1 


gain% ( 12 ) 

= 9 

PP ( 12 ) 

= 2: offset(12) 

= -1 


gain% (13 ) 

= 9 

PP (13) 

= 2: offset (13 ) 

= -1 


gain% (14 ) 

= 9 

PP (14) 

= 2: offset (14) 

= -1 


gain% (15 ) 

= 9 

PP (15) 

= 2: offset (15 ) 

= -1 



' A/D conversion number 
' A/D conversion start channel 
' A/D conversion stop channel 


26 



FUN% = 3 ' FUNCTION 3 

CALL PCL818HG (FUN% , SEG param% ( 0 ) ) ' Func 3 : Hardware initialization 

IF param% ( 45 ) <> 0 THEN PRINT "DRIVER INITIALIZATION FAILED !": STOP 

FUN% = 4 ' FUNCTION 100 

CALL PCL818HG (FUN% , SEG param%(0)) ' Func 100 : A/D initialization 

IF param% ( 45 ) <> 0 THEN PRINT "A/D INITIALIZATION FAILED STOP 

FUN% = 12 

CALL PCL818HG(FUN% , SEG param%(0)) ' Func 100 : D/A initialization 

IF param% (45) <> 0 THEN PRINT "D/A INITIALIZATION FAILED !": STOP 


END SUB 

SUB crunch (P(), 0(), probe%) 




REM 

★ ★ 

7hole probe decomposition algorithm 


★ ★ 

REM 

* * 

Input : p() array... 


★ ★ 

REM 

•k * 

P ( 1 ) . . P ( 7 ) are the pressures from the 7 

hole probe 

+ ★ 

REM 

* * 

P(8) is a normalizing pressure., tunnel 

Q 

* * 

REM 

★ * 

Output : o() array... 


★ * 

REM 

k ★ 

0(1) = Pstatic 


* * 

REM 

* * 

0(2) = Ptotal 


* * 

REM 

★ -k 

0(3) = alpha (degrees) 


+ * 

REM 

k k 

0(4) = beta (degrees) 


* * 




DIM NDpres ( 8 ) , gc ( 2 ) , pxt ( 2 ) 

REM Initialize 

alpha = -999 
beta = -999 

REM Normalize pressures . 

FOR j% = 1 TO 8 

NDpres (j%) = P ( j % ) / P(8) 

NEXT j% 

REM Compute the coefficients. 

PMAX = max (NDpres ( 1 ) , NDpres (2), NDpres (3), NDpres (4), NDpres (5), NDpres (6), 
NDpres (7) ) 

REM Default situation is zone 7 
ndhol% = 7 

PBAR = (NDpres (1) + NDpres (2) + NDpres (3) + NDpres (4) + NDpres (5) + NDpres (6)) / 

6 ! 

PDIFF = NDpres (7) - PBAR 
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CA1 = (NDpres (4) - NDpres(l)) / PDIFF 

CA2 = (NDpres (3) - NDpres (6)) / PDIFF 

CA3 = (NDpres (2 ) - NDpres(5)) / PDIFF 

cpbeta =1! / SQR(3!) * (CA2 + CA3) 
cpAlpha = CA1 + (CA2 - CA3 ) / 2! 


REM Check to see if flow angularity is outside zone 7 calibration range 

IF cpAlpha > 4.5 OR cpAlpha < -4.5 OR cpbeta > 4.5 OR cpbeta < -4.5 THEN 21 

REM look up alpha and beta from table 

ix = cpAlpha *2+10: jx = cpbeta *2+10 

iy% = INT(ix): jy% = INT(jx): iz = ix - iy%: jz = jx - jy% 

alpha = tablela (probe% , iy% + 1 , jy% +1) *iz*jz+ tablela (probe% , iy%* 3Y% + 

1) * (1 - iz) *jz+ tablela (probe% , iy%, jy%) * (1-iz) * (1-jz) + 

tablela (probe% , iy% + 1 , jy%) * iz * (1-jz) 

beta = table2a(probe%, iy% + 1, jy% + 1) * iz * jz + table2a (probe% , iy% , jy% + D 
* (1 - iz) * jz + table2a (probe% , iy%, jy%) * {1 - iz) * (1 - jz) + table2a (probe% , 
iy% + 1, jy%) * iz * (1-jz) 

GOTO 25 

REM sectoring algorithm - use separated flow coefficients 
21 : 

REM Zone 1 . 

IF (NDpres (1) = PMAX) THEN 
ndhol% = 1 

PDIFF = NDpres (1) - (NDpres (2) + NDpres (6)) / 2! 

cpAlpha = (NDpres (1) - NDpres (7)) / PDIFF 
cpbeta = (NDpres (6) - NDpres (2)) / PDIFF 
PBAR = (NDpres (2) + NDpres (6)) / 2! 

REM Zone 2 . 

ELSEIF (NDpres (2) = PMAX) THEN 
ndhol% = 2 

PDIFF = NDpres (2) - (NDpres (1) + NDpres (3)) / 2! 
cpAlpha = (NDpres (2) - NDpres (7)) / PDIFF 

cpbeta = (NDpres (1) - NDpres (3)) / PDIFF 
PBAR = (NDpres (1) + NDpres (3)) / 2! 

REM Zone 3 . 

ELSEIF (NDpres (3) = PMAX) THEN 
ndhol% = 3 

PDIFF = NDpres (3) - (NDpres (2) + NDpres (4)) / 2! 

cpAlpha = (NDpres (3) - NDpres (7)) / PDIFF 
cpbeta = (NDpres (2) - NDpres (4)) / PDIFF 
PBAR = (NDpres (2) + NDpres (4)) / 2! 

REM Zone 4 . 

ELSEIF (NDpres (4) = PMAX) THEN 
ndhol% = 4 

PDIFF = NDpres (4) - (NDpres (3) + NDpres{5)) / 21 
cpAlpha = (NDpres (4) - NDpres (7)) / PDIFF 
cpbeta = (NDpres (3) - NDpres (5)) / PDIFF 
PBAR = (NDpres (3) + NDpres (5)) / 2! 
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REM Zone 5 . 

ELSEIF (NDpres ( 5 ) = PMAX) THEN 
ndhol% = 5 

PDIFF = NDpres (5) - {NDpres (4) + NDpres (6)) / 2! 
cpAlpha = (NDpres (5) - NDpres (7)) / PDIFF 
cpbeta = (NDpres (4) - NDpres (6)) / PDIFF 
PBAR = (NDpres (4) + NDpres (6)) / 21 
REM Zone 6 . 

ELSEIF (NDpres (6) = PMAX) THEN 
ndhol% = 6 

PDIFF = NDpres (6) - (NDpres (5) + NDpres (1)) / 21 
cpAlpha = (NDpres ( 6 ) - NDpres(7)) / PDIFF 
cpbeta = (NDpres (5) - NDpres (1)) / PDIFF 
PBAR = (NDpres (5) + NDpres (1)) / 2! 

REM Zone 7 . 

ELSE 
GOTO 25 
END IF 

REM inside calibration space? 

’PRINT "SECTOR ndhol% / ’’CpAlpha cpAlpha, "CpBeta =” ; cpbeta; 

* LINE INPUT a$ 

IF cpAlpha > 2 OR cpAlpha < 0 OR cpbeta > 2 OR cpbeta < -2 THEN 25 


REM look up alpha and beta from table - sectored coefficients 
ix = cpAlpha *4+1: jx = cpbeta *2+5 

iy% = INT(ix): jy% = INT(jx): iz = ix - iy% : jz = jx - jy% 

alpha = tablelb <probe% , ndhol%, iy% + 1 # jy% + 1) *iz*jz + tablelb (probe% , 

ndhol%, iy% # jy% + 1) * (1 - iz) * jz + tablelb (probe% f ndhol%, iy%, jy%) * (1 - i 2 ) 

* (1 - jz) + tablelb (probe% , ndhol%, iy% + 1, jy%) * iz * (1 - jz) 

beta = table2b(probe%, ndhol%, iy% + 1, jy% +1) * iz * jz + table2b (probe% , 

ndhol% / iy%, jy% +1) * (1-iz) *jz+ table2b (probe% , ndhol%, iy%, jy%) * (1 - iz) 

* (1 - jz) + table2b(probe%, ndhol%, iy% + 1, jy%) * iz * (1 - jz) 

REM Have we determined flow angularity? If not error! 

IF ABS (alpha) > 90 OR ABS(beta) > 90 THEN 260 


25 : 


REM Null resiudal flow angularity 

IF probe% = 0 THEN alpha = alpha: beta = beta - 3 
IF probe% = 1 THEN alpha = alpha: beta = beta - 3 
IF probe% = 2 THEN alpha = alpha: beta = beta - 2 


REM Determine Static and Total Pressure Coefficients 
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IF alpha >45 THEN alpha =45 
IF alpha < -45 THEN alpha = -45 
IF beta > 45 THEN beta = 45 
IF beta < -45 THEN beta = -45 


ix = 10 + alpha / 5: jx = 10 + beta / 5 

iy% = INT(ix): iz = ix - iy% : jy% = INT(jx): jz = jx jy% 


cptot 
ndhol% - 

(1 - iz) 
cpsta 
ndhol% - 
(1 - iz) 


table3 (probe%, ndhol% - 1, iy% + 1- jy% + D * iz * i z + table3 (probe% 
1, iy% , jy% +1) * (1 - iz) *jz+ table3 (probe%, ndhol% - 1, iy%- iy%) 
*'(1 - jz) + table3 (probe% , ndhol% - 1, iy% + 1. DY%) * iz * (1 “ ^ z) 

: table4 (probe% , ndhol% - 1, iy% + 1- jy% + 1) * iz * jz + table4 <probe% 
1. iy%, jy% + 1) * (1 - iz) *jz + table4 (probe%, ndhol% - 1. iy%- :y%) 
* (1 - jz) + table4 (probe% , ndhol% - 1, iy% + 1» jy^) * iz * d _ j z ) 


PTOT = NDpres (ndhol% ) * (1! - cptot) + PBAR * cptot 

PSTA = PBAR * (1 + cpsta) - NDpres (ndhol% ) * cpsta 


tot = PTOT * P (8) 
sta = PSTA * P ( 8 ) 

Q = tot - sta 

REM * make static presure correction 


qc { 0 ) = .9 
qc ( 1 ) = .9 
qc ( 2 ) = .9 
pxt(0) =1.09 
pxt(l) = 1.11 
pxt(2) = 1.12 


tot = tot * pxt (probe% ) 

Q = Q / qc(probe%) 

sta = tot - Q 

0(1) = sta 
0(2) = tot 
0(3) = alpha 
0(4) = beta 

GOTO 270 

REM Notify user of errors. 

260 PRINT "ERROR : Can't resolve flow angularity" 
BEEP 

0(1) = -9.9 
0(2) = -9.9 
0(3) = -9.9 
0(4) = -9.9 
270 : 
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END SUB 


SUB DebugProbes 

SCREEN 0: WIDTH 80, 43 

CLS 

PRINT "PROBE DEBUG PRESS ANY KEY TO STOP" 

PRINT 

PRINT " PROBE PRESSURE" 

PRINT " 

PRINT 

VIEW PRINT 5 TO 43 
Talk8400 ("ODO") 


1101 : 

SOUND 1000, 1 
Talk.8400 ("ADI 11") 
tmx = TIMER + 1 

11102 : IF TIMER < tmx THEN 11102 


Talk8400 ( ” OD1 1") 

LOCATE 5 , 1 
FOR i% = 1 TO 26 

PRINT USING " ### +###.## psf "; i%; pktXdat(i%) * 144; 

n% = pktXdat(i%) * 144 / 2 
IF n% > 40 THEN n% = 40 
IF n% < -40 THEN n% = 40 

IF n% > 0 THEN PRINT STRING$(n%, "#"); TAB (78); " " 

IF n% <= 0 THEN PRINT STRINGS ( -n% , "x”); TAB (78); " ” 

IF i% = 1 OR i% = 2 OR i% = 9 OR i% = 16 OR i% = 23 THEN PRINT STRINGS (70, "-") 
NEXT i% 


IF INKEYS = * " THEN 1101 
CALL ibonl(psi%, 0) : psi% = 0 

LINE INPUT "PRESS ENTER TO RETURN TO MAIN MENU"; a$ 

VIEW PRINT 1 TO 43 
SCREEN 0: WIDTH 80, 43 

CLS 

END SUB 
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SUB DecodeBuffer (buffer$) 


Pktlnfo - Packet Header Information array - contains 

i) 

rc 

= 

Response Code 

2) 

rt 

= 

Response Type 

3) 

msglen 

= 

Message Length 

4) 

RetValue ! 

= Returned Value OR Error/Conf. Code 

5) 

NRows 

— 

number of rows 

6) 

NCols 

= 

number of columns 

7) 

Msno 

— 

Measurement Set Number 

8) 

Nvals 

= 

Number of values 

9) 

CRS 

= 

Cluster, Rack, Slot address 

10) 

Utype 

= 

Unit type 

11) 

Tblno 


Table Number 

12) 

Nf rames 

= 

Number of Frames per Data Point 

13) 

yr 

= 

year 

14) 

mo 

= 

month 

15) 

day 

= 

day 

16) 

hr 

= 

hour 

17) 

mn 

= 

minute 

18) 

sec 

= 

second 

19) 

msec 

= 

milli-seconds 


IF LEN (buf f er$ ) < 23 THEN PRINT "GARBLED BUFFER ERROR": GOTO exitsub 


PKTfinished = 0 

FOR i% = 1 TO 19: Pktlnfo (i%) = 0: NEXT i% 
pktmsg$ = ■ " 

■ Get current response code, response type, and packet length, 
rc = ASC (MID$ (buf f er$ , 1, 1)): Pktlnfo(l) = rc 
rt = ASC (MID$ ( buf f er$ , 2, 1)): PktInfo<2) = rt 

msglen = ASC (MID$ (buf f er$ , 3, 1)) * 256! + ASC (MID$ (buff er$ , 4, 1)): 
msglen ' msg length 

SELECT CASE rt 
CASE IS = 4 

GOSUB ACKpacket 
CASE IS = 9 

GOSUB SinglelEEEf loatPacket 
CASE IS = 19 

GOSUB StreamlEEEf loatPacket 
CASE IS = 33 

GOSUB ArraylEEEf loatPacket 
CASE IS = 35 
CASE IS = 128 


Pktlnfo ( 3 ) 
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GOSUB errorPacket 
CASE ELSE 

PRINT : PRINT : PRINT ■ Unknnown Packet type recieved : 
PRINT : LINE INPUT " Press ENTER to continue"; n$ 

END SELECT 


GOTO exit sub 


' Error Packet: rt = 128, and that indicates that 

some sort of error occured. See PktMsg$ for the 
error message returned by the system. 


errorPacket : 

GOSUB ACKpacket 

PRINT : PRINT " ERROR : ” ; pktmsg$ : PRINT : PKTf inished 
RETURN 


' Acknowledgement Packet: rt = 4 everything O.K. 


ACKpacket : 

ErConfCode& = CVL (MID$ (buffer $ , 8, 1) + MID$ (buf fer$ 7 
MID$ (buffer$, 5, 1)) 

RetValue! = ErConfCode& 

Pktlnfo (4) = ErConfCode& 
pktmsg$ = MID$ (buffer$, 9) 

PKTf inished = -1 
RETURN 


Single IEEE Floating Point Number Packet: rt = 9 


SinglelEEEf loatPacket : 


PktInfo(4) = CVS (MID$ (buf fer$ , 8, 1) + MID$ (buf f er$ , 7 1) 

MID$ (buf fer$ , 5, 1)) 
pktmsg$ = MID$ (buffer$, 9) 

PKTf inished = -1 
RETURN 


Stream Data, IEEE Single Precision: rt = 19 


type = “ ; rt 


-1 


1) + MID$ (buffer$, 6, 1) + 


+ MID$ (buffer$, 6, 1) + 


StreamlEEEf loatPacket : 
GOSUB DecodeStreamHeader 
j% = 25 
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FOR i% = 1 TO Nvals . . n 

pktXdat (i%) = CVS (MID$ (buf fer$ , j% + 3, 1) + MID$ (buf fer$ , 3% + 2, 
MID$ (buf fer$ , j% + 1, 1) + MID$ (buf fer$ , j%, 1) ) 

j% = j% + 4 
NEXT i% 


1 ) + 


PKTfinished = -1 
RETURN 


1 Array Data, IEEE Float : rt - 33 


ArraylEEEfloat Packet : ' 


Nrows = ASC (MID$ (buf f er$ , 5, 1) ) 
Ncols = ASC (MID$ (buf f er$ , 7, 1) ) 


* 256! + ASC (HID$ (buffer$, 6, 1)) 

* 256! + ASC (MID$ (buffers, 8, 1) ) 


Pktlnfo (5) = Nrows 
Pktlnfo (6) = Ncols 


index% = 1 
k% = 9 


FOR i% = 1 TO Nrows 
FOR j% = 1 TO Ncols 
pktXdat ( index% ) = 
MID$ (buffers, k% + 1 , 
k% = k% + 4 
index% = index% + 
NEXT j% 

NEXT i% 


CVS (MID$ (buffers, k% + 3, 1) 
1) + MID$ (buf fer$ , k%, 1)) 


+ MID$ (buffers, k% 


2 , 1 ) + 


PKTfinished = -1 
RETURN 






’ decodearrayheader : *************** 


DecodeArrayHeader : 


Nrows 

Ncols 


ASC (MID$ (buffers, 5, U) * 256! + ASC (MID$ (buf f er$ , 6, 1)) 
ASC (MID$ (buffers, 7, 1)) * 256! + ASC (MID$ (buffers, 8, 1) ) 


Pktlnfo (5) = Nrows 
Pktlnfo (6) = Ncols 

pktmsgS = MID$ (buf fer$ , 9, 40) 

RETURN 
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decodes treamheader : 


DecodeStreamHeader : 

we have already decoded response code, response type, 
' and response length 


Msno = ASC (MID$ (buf fer$ , 5, 1)) * 256! + ASC (MID$ (buf f er$ , 6, 1)) 
Nvals = ASC (MID$ (buffer$, 7, 1)) * 256! + ASC (MID$ (buf fer$ , 8, 1)) 
crs = ASC (MID$ (buffers, 9, 1)) * 100 + ASC (MID$ (buf fer$ , 10, 1)) * 10 
crs = crs + ASC (MID$ (buf fer$ , 11, 1)) 

Utype = ASC (MID$ (buffers, 12, 1)) 

Tblno = ASC (MID$ (buffers, 13, 1)) 

Nframes = ASC (MID$ (buf f er$ , 14, 1)) 

yr = ASC (MID$ (buffers, 15, 1)) 
mo = ASC (MID$( buffers, 16, 1)) 
dy = ASC (MID$ (buffers, 17, 1)) 
hr = ASC (MID$( buffers, 18, 1)) 

mn = ASC (MID$( buffers, 19, 1) ) 
sc = ASC (MID$ (buffers, 20, 1)) 


msec - ASC (MID$ (buffers, 21, 1)) * 256! + ASC (MID$ (buf fer$ , 22, 1) ) 


Pktlnfo {7 ) 
PktInfo(8) 
Pktlnfo (9) 
Pktlnfo ( 10 ) 
Pktlnfo ( 11 ) 
Pktlnfo (12) 
Pktlnfo (13) 
Pktlnfo ( 14 ) 
Pktlnfo ( 15 ) 
Pktlnfo ( 16 ) 
Pktlnfo (17) 
Pktlnfo { 18 ) 
Pktlnfo ( 19 ) 


= Msno 
= Nvals 
- crs 
= Utype 
= Tblno 
= Nframes 
= yr 
= mo 
= dy 
= hr 
= mn 
= sc 
= msec 


RETURN 


exitsub: 


END SUB 


SUB GetPosition 
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I - 

'Parallel Read Encoder Position for MetraByte PIO-12 and AR Absolute 
'encoders. The ACT array is used to determine which encoders will be 
'Read. The output will be in the arrays ERR and PLC. The method 
•used to reduce the parallel data will be determined by the status of 
'the SP parameter for each encoder (either Binary or BCD). 

•PARALLEL READ ENCODER POSITION 

DIM pic (4) 

' Set-up ' s for MetraByte Parallel Ports 


’Converse in parallel 

OUT CTRL, CONTROL 1 Set up MetraByte 

ERP = false 

FOR R% = 0 TO 4 
IF ACT ( R% ) <> 0 THEN 

DEV = 2 A (R% + 3) ‘Select which encoder to poll for position and error data 

OUT PC, 4 + DEV + 3 'Strobe High (4), AO and Al High (3), Device Selected 
OUT PC, DEV + 3 'Strobe Low, etc. 

FOR i% = 0 TO delay: NEXT i% 

LSB = INP(PA) 'Get LSB when AO and Al high 
OUT PC , DEV + 1 ' Drop Al Low 

FOR i% = 0 TO delay: NEXT i% 

LSB2 = INP(PA) 'Get LSB2 when AO high and Al low 
OUT PC, DEV 'Drop AO Low as well 

FOR i% = 0 TO delay: NEXT i% 

MSB2 = INP(PA) 'Get MSB2 when AO and Al low 
OUT PC, DEV + 2‘ Raise Al high 

FOR i% = 0 TO delay: NEXT i% 

MSB = INP(PA) 'Get MSB when AO low and Al high 
OUT PC, 4 + DEV + 3 'Bring up all lines 

'Convert data to position and error information 
IF ESB$(R%, 6) = "0" THEN 
'Data interpreted as Binary 

SIGN = 1: IF (LSB AND 16) = 16 THEN LSB = LSB - 16: SIGN - -1 

IF (LSB AND 64) = 64 THEN ERP = true: LSB = LSB - 64 

IF (LSB AND 128) = 128 THEN ERP = true: LSB = LSB - 128 

plc(R%) = ((LSB / 65536!) + (LSB2 / 512) + MSB2 + (MSB * 256)) * SIGN 
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ELSE 


Data, interpreted as Binary Coded Decimal 

SIGN =1: IF {MSB AND 16) = 16 THEN MSB = MSB - 16: SIGN = -1 
IF {MSB AND 64) =64 THEN ERP = true: MSB = MSB - 64 
IF {MSB AND 128) = 128 THEN ERP = true: MSB = MSB - 128 
pic (R% ) = VAL (HEX$ (MSB) ) * 100 + VAL {HEX$ (MSB2 ) ) 

pic (R% ) = pic (R% ) + VAL (HEX$ (LSB2 ) ) * .01 + VAL ( HEX$ ( LSB ) ) * .0001 
pic (R% ) = pic ( R% ) * SIGN 

END IF 

END IF 

IF ABS (pic (R% ) ) > 69.999 THEN ERP = true: plc(R%) = -999 
NEXT R% 

IF ERP = true THEN SOUND 5000, 1 

IF Pled ) = -999 THEN ELSE tX = plc(l) 'X (longitudinal) 

IF plc(0) = -999 THEN ELSE ty = plc{0)' Y (horizontal) 

IF pic (2 ) = -999 THEN ELSE Tz = plc(2) ' Z (vertical) 

cx = POS(O): cy = CSRLIN: VIEW PRINT 1 TO 43 

COLOR 14, 1: LOCATE 3, 1: PRINT USING "TRAVERSE POSITION : HORZ +###.###in 
+###.### in LONG +###.### in " ; ty; Tz; tX; TAB (80); 

VIEW PRINT 4 TO 41: LOCATE cy, cx 
COLOR 7,0 


END SUB 

SUB Init 

SCREEN 0 
WIDTH 80, 43 

PRINT 

PRINT 

PRINT "DONE.” 

PRINT 

PRINT M SLOWSURVEY . BAS - WAKE SURVEY MODULE - 1995 - T . Takahashi ” 
PRINT 

PRINT “ Initializing....” 

CHDIR "c : \qbasic\pcdas 11 

OPEN "I”, 1, N pc-das . ini " 

LINE INPUT #1, pt$ 

LINE INPUT #1, tes t$ 

CLOSE 1 


VERT 
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IF RIGHT$(pt$, 1) 

1 ) 


\" AND RIGHT$ (pt$ , 2 ) <> 


:\" THEN pt$ = LEFT$(pt$, LEN(pt$) 


CHDIR pt$ 

IF RIGHT$ (pt$ , 1) <> "V THEN pt$ = pt$ + "V 
PRINT " Reading 7HP Calibration Data" 


END SUB 

FUNCTION keystroke% (0$) 


1000 : 

CALL GetPosition 

a$ = INKEY $ 

IF a$ = "" THEN 1000 


COLOR 7 , 0 : PRINT a$ 
k% = INSTR (0$ , a$) 

IF k % = 0 THEN BEEP: COLOR 4, 0: PRINT "KEYSTROKE ERROR! - VALID COMMANDS 

COLOR 7 1 

keystroke% = k% 

END FUNCTION 


SUB Main 


1 : 

IF motor = 1 THEN PRINT " MOTOR CONTROL DISABLED" 

Menu ("SURVEY ALIGN PROBE_DEBUG 
COLOR 15 

PRINT "WAKE SURVEY CONTROL" 

COLOR 7 

PRINT " ENTER COMMAND : " ; 

op = keystroke ( "SsAaPpEe" ) 

SELECT CASE op 

CASE 1, 2: CALL SloSurvey 
CASE 3, 4: CALL PositionMenu 
CASE 5, 6: CALL DebugProbes 

CASE 7, 8: GOTO done 

END SELECT 
GOTO 1 


" ; 0 $ : 


EXIT” ) 
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done: 


END SUB 

FUNCTION max (a, b, c, d, e, f, g) 
n = 0 

IF a > b THEN n = a ELSE n = b 

IF c > n THEN n = c 

IF d > n THEN n = d 

IF e > n THEN n = e 

IF f > n THEN n = f 

IF g > n THEN n = g 

max = n 

END FUNCTION 

SUB Menu (m$) 

cx = POS(O): cy = CSRLIN 

VIEW PRINT 1 TO 43 

LOCATE 1, 1: COLOR 7, 1: PRINT TAB (37); "PC-DAS"; TAB (80); 

LOCATE 2,1: 

flag% = 32 

FOR i% = 1 TO LEN ( m$ ) 

IF f lag% = 32 THEN COLOR 15, 7 ELSE COLOR 0, 7 
f lag% = ASC (MID$ (m$ , i%, 1)) 

PRINT CHR$ ( f lag% ) ; 

NEXT i% 

PRINT TAB (80) ; 

LOCATE 42, 1: COLOR 0, 3 

PRINT "PATH: pt$; TAB(30); "TEST: test$; TAB(60); TIME$ ; " DATE$ ; TAB(79); 


VIEW PRINT 4 TO 41 
IF cy < 4 THEN cy = 4 
LOCATE cy, cx 
COLOR 7, 0 

END SUB 


FUNCTION min (a, 
n = 0 

IF a < b THEN n 
IF c < n THEN n 
IF d < n THEN n 
IF e < n THEN n 
IF f < n THEN n 
IF g < n THEN n 


b/ c, d, e, f, g) 


a ELSE n = b 

c 

d 

e 

f 

g 


min = n 
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END FUNCTION 


SUB moveto (Hnew, Vnew, Lnew) 


'Variable Reference 
•HPos,VPos, tX 
■HOld, VOld, LOld 
' HNew , VNew , LNew 
■ HHome , VHome , LHome 
' HDist, VDist, LDist 
’ HStep , VStep , LS tep 


Current Position (Horizontal, Vertical, and Longitudinal) 
= Last (Previous) Position 
= New (Given/Expected) Position 
= Home Position 

= Distance from Current to New (ex. HD=HN-Hor zC 
= Compumotor Steps to achieve Distance 


PRINT USING 


Move to : +##.### +##■### +##-###"; Hnew; Vnew; Lnew 


IF motor% = 1 THEN GOTO 12345 
■ make sure the traverse has stopped 


tO = TIMER 

' Where are we? (we already know from the call to still moving) 
loopback: 

1 How far do we need to go? 

HDist = Hnew - ty 
VDist = Vnew - Tz 
LDist = Lnew - tX 

'Compute Motor Steps 
Hstep& = -HDist * Hscale 
Vstep& = VDist * VScale 
Lstep& = LDist * LScale 

' PRINT Hstepfc ; Vstep&; Lstep&c, ; 

’Which Way are we going? 
vl = SGN (VDist) 
v2 = SGN (HDist) 

' move specified distance minus last little bit 
IF ABS (Vstep&) >= 199 THEN Vstep& = Vstep& - vl * 30 
IF ABS (Hstepi) >= 1299 THEN Hstep& = Hstepfc + v2 * 10 

Vdir = vl 

CALL TalkCompumotor (Hstepfc, Vstepk, Lstep&) 

IF ABS(Hstep&) <= 10 AND ABS(Vstep&) <= 40 THEN GOTO 12345 


40 



REM * Wait one second for the compumotor action to begin 


tO = TIMER + 1 
DO 

LOOP UNTIL (TIMER > tO) 
CALL StillMoving 
GOTO loopback 
12345 : 

BEEP 
END SUB 


SUB PositionMenu 
2 : 


Menu ("MOVE_TO UP DOWN LEFT RIGHT EXIT") 

COLOR 15 

PRINT " ALIGN SURVEY RIG" 

COLOR 7 

PRINT " ENTER COMMAND : " ; 

op = keystroke ( "MmUuDdLlRrEe" ) 

SELECT CASE op 
CASE 1, 2: GOTO SmartMove 
CASE 3, 4: GOTO Up 
CASE 5, 6: GOTO Down 
CASE 7, 8: GOTO Left 
CASE 9, 10: GOTO Right 
CASE 11, 12: GOTO 3 

END SELECT 

GOTO 2 

H = ty: V = Tz: 1 = tX 


Up: 

IF Tz < VHiLim THEN CALL moveto(ty, Tz + 1, tX) 
GOTO 2 

Down: 

IF Tz > VLoLim THEN CALL moveto(ty, Tz - 1, tX) 
GOTO 2 

Left : 

IF ty > HLoLim THEN CALL moveto(ty + 1, Tz, tX) 
GOTO 2 

Right : 
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IF ty < HHiLim THEN CALL moveto(ty - 1, Tz, tX) 

GOTO 2 

SmartMove : 

INPUT " Enter New Position (H,V,L) : " ; H, V, 1 

IF H < HLoLim OR H > HHiLim OR V < VLoLim OR V > VHiLim OR I < LLoLim OR 1 
THEN 


BEEP 

COLOR 4 , 0 






PRINT “ 

ERROR 

, ! KEEP 

RIG WITHIN 

LIMITS ! " 

PRINT " 

HORZ 

from " ; 

HLoLim; “ 

to " ; 

HHiLim 

PRINT " 

VERT 

from " ; 

VLoLim; ■ 

to " ? 

VHiLim 

PRINT " 
PRINT 
GOTO 2 
END IF 

LONG 

from " ; 

LLoLim; " 

to ” ; 

LHiLim 


PRINT " Moving Traverse (press any key to interrupt) " 

CALL move to (H, V, 1) 

BEEP 

GOTO 2 


3 : 

END SUB 

SUB readanalog 

DIM param%(60), DAT%(32000), gain%(15), nn{15) , pp(15), offset(15) 


param% ( 0 ) = 

0 

' Board number 

param% ( 1 ) = 

&H300 

' Base I/O address 

par am % (4 ) = 

2 

' IRQ level : IRQ2 

param% { 6 ) = 

20: param%(5) = 20: 

xmax = 2000 ' sample at 500hz 

param% ( 7 ) = 

0 

* Trigger mode, 0 : pacer trigger 

param% ( 8 ) = 

0 

* Non-cyclic 

param% ( 10 ) 

= VARPTR ( DAT% ( 0 ) ) 

1 Offset of A/D data buffer A 

param% (11) 

= VARS EG ( DAT% ( 0 ) ) 

1 Segment of A/D data buffer A 

param% ( 12 ) 

= 0 

' Data buffer B address, if not used, 

param% { 13 ) 

= 0 

' must set to 0. 

param% ( 14 ) 

= 4096 

' A/D conversion number 

param% (15) 

= 0 

' A/D conversion start channel 

param% ( 16 ) 

i 

= 15 

' A/D conversion stop channel 

param% ( 17 ) 
param% ( 18 ) 

= &HFF ' gain table 
= VARPTR (gain% (0) ) 

1 Offset of gain table buffer 


> LHiLim 
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param% (19) 

= 

VARSEG(gain% (0) ) ' Segment of gain table buffer 

gain% ( 0 ) 

= 

0: 

pp ( 0 ) = 

10: 

offset(0) = 

-5: ' +/- 5v external input 

gain% (1) 

= 

0: 

pp ( 1 ) = 

10: 

offset(l) = 

-5: 1 Pa 

gain% (2 ) 

= 

8: 

PP ( 2 ) = 

20: 

offset(2) = 

-10: ' Pr-Pa 

gain% (3 ) 

= 

8: 

PP ( 3 ) = 

20 : 

offset(3) = 

-10: 1 Pr-Ps 

gain% (4 ) 

= 

4: 

pp ( 4 ) = 

10: 

offset(4) = 

0 : ' turntable 

gain% (5) 

= 

0: 

PP ( 5 ) = 

10: 

offset(5) = 

-5: 1 inclinometer 

gain% ( 6 ) 

= 

6: 

PP ( 6 ) = 

.1: 

offset(6) = 

0 : ’ PTplenuin barocel 

gain% (7 ) 

= 

6: 

PP ( 7 ) = 

.1: 

offset(7) = 

0: ' Tplenum thermocouple 

gain% (8) 

= 

4: 

pp ( 8 ) = 

10: 

offset(8) = 

0 : ’ CFM plenum 

gain% (9 ) 

= 

9: 

PP ( 9 ) = 

2 : 

offset(9) = - 

-1: ' PS forebody 

gain% (10) 

= 

: 9 

pp(10) 

= 2 

: offset(10) 

= -1 

gain% (11) 

= 

9 

PP(ll) 

= 2 

: offset ( 11 ) 

= -1 

gain% ( 12 ) 

= 

9 

PP (12) 

= 2 

: offset ( 12 ) 

= -1 

gain% (13) 

- 

9 

PP (13) 

= 2 

: offset(13) 

= -1 

gain% (14 ) 

= 

9 

PP (14) 

= 2 

: offset(14) 

= -1 

gain% (15) 

= 

9 

PP (15) 

= 2 

: offset(15) 

= -1 


FUN% = 5 1 FUNCTION 5 

CALL PCL818HG ( FUN% , SEG parain%(0) ) ' Func 5 : Pacer trigger A/D 

IF param% ( 45 ) <> 0 THEN PRINT "A/D INTERRUPT DATA TRANSFER FAILED !": 

FOR i% = 0 TO 15 
volts (i%) = 0: nn(i%) = 0 
NEXT i% 

FOR i% = 0 TO param% ( 14 ) - 1 
ix% = i% MOD 16 

volts (ix%) = volts (ix%) + (pp(ix%) * DAT% ( i% ) / 4096! + offset (ix%) ) 
nn ( ix% ) = nn ( ix% ) + 1 
NEXT i% 

FOR i% = 0 TO 15 
volts(i%) = volts(i%) / nn(i%) 

NEXT i% 

END SUB 

SUB ReadCal 


PRINT "READING CALIBRATION MATRICIES " 
OPEN "I", 1, ‘'c:\qbasic\pcdas\shpcal.dat" 
FOR probe% = 0 TO 2 
PRINT " PROBE " ; probe% 

FOR ix% = 1 TO 20 
FOR jx% = 1 TO 20 
FOR ndhol% = 0 TO 6 
INPUT #1, a, b, c, d, e, f 
table3 (probe% , ndhol%, ix%, jx%) = e 
table4 (probe% , ndhol%, ix% , jx%) = f 


STOP 
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NEXT ndhol% 

NEXT jx% 

NEXT ix% 

ndhol% = 7 

FOR cpAlpha = -4.5 TO 4.5 STEP .5 
FOR cpbeta = -4.5 TO 4.5 STEP .5 
INPUT #1, a, b, c, d, e, f 
x% = cpAlpha *2+10 
y% = cpbeta *2+10 
tablela(probe%, x%, y%) = e 
table2a (probe% , x%, y%) = f 
NEXT cpbeta 
NEXT cpAlpha 

FOR ndhol% = 6 TO 1 STEP -1 
FOR cpAlpha = 0 TO 2 STEP . 2 5 
FOR cpbeta = -2 TO 2 STEP .5 
INPUT #1, a, b, c, d, e, f 
x% = cpAlpha *4+1 
y% = cpbeta *2+5 

tablelb (probe% , ndhol%, x%, y%) = e 
table2b(probe%, ndhol%, x%, y%) = f 
NEXT cpbeta 
NEXT cpAlpha 
NEXT ndhol% 

NEXT probe% 

CLOSE 1 

END SUB 

SUB readscales 
DIM s (10) 

REM * PIO 96 base address 
basead2 = &H2D0 

REM * reset 
resetmux = &H4 
inhibitl = &H2 
inhibit2 = &H1 

REM * set up PIO 96 to talk to digital MUX 
OUT basead2 + 3 , &H9B 
OUT basead2 + 7, &H90 

OUT basead2 +5,0 ' reset scales 

REM * reset mux 
OUT basead2 + 5, resetmux 
FOR i = 0 TO 9 9 : NEXT 
OUT basead2 +5, 0 
REM * read data 
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FOR n% = 0 TO 9 
OUT basead2 + 5, inhibit2 
OUT basead2 + 5, 0 
b4 = INP(basead2 + 4) 
b3 = INP (basead2 + 2) 
b2 = INP (basead2 + 1) 
bl = INP (basead2 ) 

OUT basead2 + 5, inhibitl 
OUT basead2 +5, 0 

IF (b4 AND 128) =128 THEN b4 = b4 - 128: SIGN = -1 ELSE SIGN = 1 
IF (bl > &H99) OR (b2 > &H99) OR (b3 > &H99) OR (b4 > &H99) THEN 
num = -999 . 99 
ELSE 

num = SIGN * ( VAL (HEX$ (b4 ) ) * 100000 + VAL ( HEX$ ( b3 ) ) * 1000 + VAL {HEX$ (b2 ) ) * 10 
+ VAL { HEX$ (bl) ) / 10) 

END IF 
s ( n% ) = num 
NEXT n% 


FOR n% = 0 TO 9 

IF s(n%) <> -999.99 THEN scale(n%) = s(n%) 
NEXT n% 


END SUB 


SUB SetUpCompumotor 

'Set-up Encoder Status Block for 5 AR W/RS422 Serial Absolute Encoders. 
’The "Status Block" is a series of array variables that contain (after 
’initialization of the encoders) the current status of each encoder. By 
’modifying these initial values, the encoders can be brought up in any 
'state. The parameter strings (P$(n)) hold the command string used to 
'communicate with the encoders. The ACT array holds either 0 or 1 for 
’each encoder. If the encoder is to be used (ie. active) the ACT varaible 
'will hold a 1 for that encoder. Otherwise ACT will hold a zero, causing 
'the initialization section to ignore that encoder. 


'ENCODER STATUS BLOCK 

'Set-up's for Parameter command strings 


P$ (0) 
P$(l) 
P$(2) 
P$(3) 
P$ ( 4 ) 
P$ (5) 
P$ (6) 
P$(7) 


"SN" ' Set 
"SE" ' Set 
"SD" ' Set 
” SM " ' Set 
"SO” 'Set 
"SF" ' Set 
”SP" ’Set 
"AR" ' Set 


encoder 

encoder 

encoder 

encoder 

encoder 

encoder 

encoder 

encoder 


Address Character (0-15) 

Error Checking (0-Off,l-On) 

Direction (0-inc. cclkwise , 1-inc . clockwise) 
Serial Duplex ( 0-Full # 1-Half >see manual<) 
Position Offset (0-cur. pos=0 . 000 , A-pos=absolute) 
Scale Factor (uffff-units digit & 4 dec. places) 
Data Format ( 0-Bin/Hex, 1-BCD/Dec [par/ser] ) 

Auto Report (0-Off,l-On) 


'Initialize each Encoder Block 


' Encoder 0 
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ESB$(0, 0) = 
ESB$(0, 1) = 
ESB$ ( 0 , 2) = 
ESB$ (0 , 3) = 
ESB$(0, 4) = 
ESB$(0, 5) = 
ESB$ ( 0 , 6) = 
ESB$ (0 , 7) = 
ESB$ ( 0 , 8) = 

' Encoder 1 
ESB$ ( 1 , 0) 
ESB$ ( 1 , 1) 

ESB$(1, 2) 
ESB$ (1 , 3) 
ESB$ ( 1 , 4) 
ESB$ ( 1 , 5 ) 
ESB$ ( 1 , 6) 

ESB$ (1 , 7) 
ESB$ ( 1 , 8) 

' Encoder 2 
ESB$ (2 , 0) 
ESB$ (2 , 1) 
ESB$ (2 , 2) 
ESB${2, 3) 
ESB$ (2 , 4) 

ESB$ ( 2 , 5) 
ESB$ (2 , 6) 

ESB$ (2 , 7) 
ESB$ (2 , 8) 

' Encoder 3 
ESB$ (3 , 0) 
ESB$ (3 , 1) 

ESB$ ( 3 , 2) 
ESB$ (3 , 3) 
ESB$ (3 , 4) 
ESB$ (3 , 5) 
ESB$ (3 , 6) 

ESB$ (3 , 7) 
ESB$ (3 , 8) 

' Encoder 4 
ESB$ (4 , 0) 
ESB$ ( 4 , 1) 
ESB$ ( 4 , 2) 
ESB$ (4 , 3) 
ESB$ (4 , 4) 
ESB$ ( 4 , 5) 
ESB$ (4 , 6) 

ESB$ (4 , 7) 
ESB$ (4 , 8) 


" 0 " ' SN Encoder 0 address is 0 

* l •' ' SE Error checking OFF 

"1" 1 SD Encoder counts increase as shaft turns clockwise 

"l" ' SM Serial HALF duplex 

"A" 'SO Position Offset OFF (ie. position is absolute) 

" 10000 ” ' SF Scale Factor is 1.0000 (ie. counts, not eng. units) 
"2“ ' SP Parallel output in BCD, Serial output in Dec. 

“1" 'AR Auto Report Mode ON 
"Horizontal “ 


" 1" ' SN Encoder 1 address is 1 

■ 1 " ' SE Error checking ON 

"1" ' SD Encoder counts increase as shaft turns clockwise 

" 1" 'SM Serial HALF duplex 

"A" 'SO Position Offset OFF (ie. position is absolute) 

“10000" ' SF Scale Factor is 1.0000 (ie. counts, not eng. units) 
" 2 " ' SP Parallel output in BCD, Serial output in Dec. 

"1" 'AR Auto Report Mode ON 
"Longitudinal “ 


"2" 1 SN Encoder 2 address is 2 

"0" 1 SE Error checking OFF 

"1" 1 SD Encoder counts increase as shaft turns clockwise 

"1" ■ SM Serial HALF duplex 

"A" ’SO Position Offset OFF (ie. position is absolute) 

”10000" 1 SF Scale Factor is 1.0000 (ie. counts, not eng. units) 
"1" 1 SP Parallel output in BCD, Serial output in Decimal 

_ . a R Auto Report Mode ON 

= "Vertical 


= “3" • SN Encoder 3 address is 3 

= " 0 " ■ SE Error checking OFF 

= "1” ' SD Encoder counts increase as shaft turns clockwise 

= "1" 'SM Serial HALF duplex 

= "A" 'SO Position Offset OFF (ie. position is absolute) 

= "10000" ' SF Scale Factor is 1.0000 (ie. counts, not eng. units) 
_ n ^ ii 1 SP Parallel output in BCD, Serial output in Decimal 
_ i ■ AR Auto Report Mode ON 
= "Encoder 3 


= ” 4 " ' SN Encoder 4 address is 4 

= " 0 “ ' SE Error checking OFF 

- n " • sd Encoder counts increase as shaft turns clockwise 

= "1" ' SM Serial HALF duplex 

= “A" 'SO Position Offset OFF (ie. position is absolute) 

= "10000" ' SF Scale Factor is 1.0000 (ie. counts, not eng. units) 

- “i" • SP Parallel output in BCD, Serial output in Decimal 

= "1" 'AR Auto Report Mode ON 

= "Encoder 4 
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'Set-up's for each encoder 


ACT ( 0 ) =1 

1 Encoder 

1 = on 

(Horizontal ) 

ACT ( 1 ) =1 

' Encoder 

1 = on 

{ Longi tudinal ) 

ACT (2) = 1 

' Encoder 

2 = on 

(Vertical ) 

ACT (3) = 0 

' Encoder 

3 = off 


ACT (4) = 0 

1 Encoder 

4 = off 



hveloc = 12 
END SUB 
SUB SloSurvey 
DIM P ( 8 ) , 0(4) 

CALL analoginit 
COLOR 15 

PRINT " WAKE SURVEY" 

COLOR 7 

PRINT “ ALIGN SURVEY RIG TO DESIRED LONGITUDINAL POSITION" 

PRINT 

PRINT ’’ OK TO CONTINUE? [Y/N] " 


op = keystroke ( "YyNn" ) 
SELECT CASE op 
CASE 1, 2: GOTO 10000 
CASE 3, 4: GOTO 20000 
END SELECT 


10000 : 

PRINT " MODIFIED FOR FLAP EDGE TEST" 

INPUT " ENTER STARTING COORDINATES ( HORZ , VERT ) : hO , vO 

INPUT " ENTER ENDING COORDINATES (HORZ, VERT) : hi, vl 

PRINT 

LINE INPUT " ENTER COMMENTS : “ ; c$ 

IF hi < hO THEN SWAP hO , hi 
IF vl < vO THEN SWAP vO , vl 


dv = vl - vO: IF dv < 2 THEN dv = 2 
dh = hi - hO: IF dh < 10 THEN dh = 10 

CALL GetPosition: 10 = tX 

PRINT " MOVING TO INITIAL SURVEY LOCATION" 

CALL moveto(h0, vO , 10) 

LINE INPUT "ENTER FILESPEC FOR WAKE SURVEY DATA f$ 

IF f$ = " ” THEN f$ = "c:\noname.dat" 
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OPEN "0", 2, f$ 


PRINT M BEGINNING SLOW SURVEY" 

Talk8400 ("SC4 0")' disable front panel echo 


WRITE #2 , test$ 

WRITE #2, DATE$ , TIME$ 


WRITE 

to 

o 

</> 




PRINT 

#2, "" 




WRITE 

#2, "X" 

h v " 

t L t 

" z " , " Q " , " TT " 

, "RH", "BARO", "Pll", 

" P15 " , 

" P16 " , 

" PI 7 " , 

" PT1 " , "PS1", 

"Al" , "B1 ” / “ P21 “ , ” P22 

“P26" , 

"P27" , 

“ PT2 " , 

" PS2 " , " A2 " , " 

B2 ” , " P31 ” , " P32 " , “ P33 

" P37 “ , 





” PT3 " , 

"PS3" , 

"A3 " , 

" B3 " 



P12 " , 

# " P23 
, "P34 


P13 " , 

, " P24 

, " P35 


vO = vO - 1 


WHILE Tz < vl 
vO = vO + 1 

REM * auto horz. domain for flap-edge test 

hveloc = 1: REM set a slow horizontal drift velocity 
hO = 0: hi = 15: dh = 15 

IF vO > 20 THEN hO = -25: hi = 0 : dh = 25 

IF vO <= 20 THEN hO = -10: hi = 5: dh = 15 

IF INT ( vO / 4) * 4 = vO THEN hO = -30: hi = 25: dh = 55 

hveloc - 10 


CALL GetPosition 

IF ABS ( ty - hO) < .1 THEN 
1 at left boundary 
CALL moveto(hO, vO , 10} 
hx = hi 
direction = 1 
ELSE 

CALL moveto(hl, vO, 10) 
hx = hO 
direction = -1 
END IF 

REM * Here is a good place to acquire non-7hp telemetry 
REM * Barocell, Tunnel Temp, Humidity 

CALL readanalog 
CALL readscales 


P14 " , 

, " P2 5 

, “ P3 6 
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qscale = scale (6) 

IF scale ( 7 ) > -110 AND scale(7) < 0 THEN tt = -scale(7) 
rh = scale(8) / 10 

pBARO = volts (1) / .1001992: IF pBARO < 13 OR pBARO > 16 THEN pBARO = 14.7 

qcse = qscale ' converted to psf 

• STATIC PLATE CORRECTION FUNCTION (after Wadcock, 7xl0ist, 1996) 

' aft static ring (#1) - nearest test section 

’ dqsp = -0.01197 - .01395 * qcse - .000145 * qcse A 2 - 3.491 E-07 * qcse 
' mid static ring (#2) 

•dqsp = .00773 - .1191 * qcse - 7.13e-5 * qcse A 2 - 3.67E-07 * qcse A 3 

• fwd static ring (#3) - nearest settling chamber 

dqsp = .00011 - .1571 * qcse - .000148 * qcse A 2 - 2.603E-07 * qcse A 3 

Qps = qcse - dqsp ' q from pitot-static on c/1 of tunnel - using static ring # 3 


HDist = hi - hO : Hstep& = -HDist * Hscale * direction: REM CALC steps 
hveloc = 1: REM return to usual horizontal speed 
IF INT(v0 / 4) * 4 = vO THEN hveloc = 2 

CALL TalkCompumotor (Hstep& , 0, 0) 

REM * wait a wee bit to let the motors spool up 
tO = TIMER + 1 
DO 

LOOP UNTIL (TIMER > tO) 

PRINT " Y Z | QPRBE | PT1 PS1 A1 B1 | PT2 PS2 A2 B2 | PT3 PS3 A3 

B3 “ 


WHILE (ABS(ty - hx) > .1) 

REM * ACQUIRE PSI-DATA HERE 
Talk8400 ("ADI 11") 


REM * Get out position while the PSI is out doing its thing 
CALL GetPosition 

REM * We can also start doing disk I/O as well 

PRINT #2, USING "+##.#!+##.#!+##.#!“; tX; ty; Tz; 

PRINT #2, USING “+##.#!+##.#!+##.###!+##•#!"; tt; * , ” ; rh; 


pBARO ; Qps; 


REM Wait 1/5 econd 
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to = TIMER + .2 
DO 

LOOP UNTIL (TIMER > tO) 


REM * get the data from the PSIs 
Talk8400 ( “OD1 " ) 

tO = Pktlnfo ( 18 ) + Pktlnfo ( 19 ) / 1000 

PRINT USING "+##.# +##.# | +### |"; ty; Tz; Qps; 

probe% = 0 : P ( 1 ) = pktXdat ( 1 ) : P ( 2 ) = pktXdat { 2 ) : P 
P ( 4 ) = pktXdat ( 4 ) : P ( 5 ) = pktXdat ( 5 ) : P ( 
P ( 7 ) = pktXdat ( 7 ) : P ( 8 ) = Qps 
CALL crunch (P(), 0(), probe%) 

PRINT USING "+## +## +## +## | “ ; 0(1) * 144; 0(2) 
PRINT #2, USING "+#.####!"; P(l); P(2); 

P ( 6 ) ; * , * ; P ( 7 ) ; 

PRINT #2, USING "+#.####!"; 0(1); 0(2); “ , " ; 

PRINT #2, USING "+##.#!"; 0(3); 0(4); 

tpO =0(1) 

probe% = 1 : P ( 1 ) = pktXdat ( 8 ) : P ( 2 ) = pktXdat ( 9 ) : P 
P ( 4 ) = pktXdat (11): P(5) = pktXdat (12): 
P(7) = pktXdat (14): P(8) = Qps 
CALL crunch ( P ( ) , O ( ) , probe% ) 

PRINT USING "+## +## +## +## | 0(1) * 144; 0(2) 

PRINT #2, USING “+#.####!“; P(l); P { 2 ) ; 

P ( 6 ) ; P ( 7 ) ; ; 

PRINT #2, USING "+#.####!"; 0(1); " , " ; 0(2); 

PRINT #2, USING "+##.#!"; 0(3); 0(4); 

tpl = 0(1) 

probe% =2: P(l) = pktXdat (15): P ( 2 ) = pktXdat (16): 
P ( 4 ) = pktXdat (18): P ( 5 ) = pktXdat (19): 

P ( 7 ) = pktXdat (21): P ( 8 ) = Qps 
CALL crunch ( P ( ) , 0 ( ) , probe% ) 

PRINT USING "+## +## +## +##"; 0(1) * 144; 0(2) * 
PRINT #2, USING "+#.####!"; P(l); " , " ; P(2); " , " ; 

P ( 6 ) ; ; P ( 7 ) ; 

PRINT #2, USING "+#.####!"; 0(1); 0(2); 

PRINT #2, USING "+##.#!"; 0(3); 0(4); M " 

tp2 =0(2) 


WEND 

BEEP 

CALL ibonl(psi%, 0) 
psi% = 0 

IF INKEY$ = CHR$ ( 27 ) THEN 30000 
WEND: REM * have finished a horizontal traverse 


30000 : 


(3) = pktXdat (3) 
6 ) = pktXdat ( 6 ) 


* 144; 0(3) ; 0(4) ; 

P ( 3 ) ; P ( 4 ) ; P(5) 


(3) = pktXdat (10) 
P(6) = pktXdat (13) 


* 144; 0(3) ; 0(4); 

P ( 3 ) ; * , ■ ; P ( 4 ) ; " , " ; P ( 5 ) 


P(3) = pktXdat (17) 
P ( 6 ) = pktXdat (20) 


144; 0(3); 0(4) 

P ( 3 ) ; " , " ; P ( 4 ) ; " , " ; P ( 5 ) 
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Talk8400 ("SC4 1")' front panel echo on 
CALL ibonl(psi% / 0): psi% = 0 


PRINT " SURVEY COMPLETE - RETURNING TO HOME POSITION” 
CALL moveto(h0, vO , 10) 


CLOSE 2 


20000 : 


hveloc = 12 
END SUB 

SUB SmartMove 

GetPosition 


INPUT "Enter New Position (H,V,L) : M ; ty, Tz , tX 

H = ty: V = Tz: 1 = tX 
CALL moveto{H, V, 1) 

BEEP 

END SUB 

SUB StillMoving 
DO 

GetPosition 

LOOP UNTIL ( ERP = false) 
checkagain: 

HOld = ty: VOld = Tz : LOld = tX 

tO = TIMER + .5: WHILE (TIMER < tO): WEND 

DO 

GetPosition 

LOOP UNTIL (ERP = false) 


IF INKEY$ <> " ” THEN GOTO KeyBreak 
IF (ty < HLoLim) THEN GOTO OutOf Bounds 

IF (ty > HHiLim) THEN GOTO OutOf Bounds 

IF (Tz < VLoLim) THEN GOTO OutOf Bounds 

IF (Tz > VHiLim) THEN GOTO OutOf Bounds 

IF (tX < LLoLirn) THEN GOTO OutOf Bounds 

IF (tX > LHiLim) THEN GOTO OutOf Bounds 
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IF ({HOld = ty) AND (VOld = Tz) AND (LOld = tX) ) THEN 
GOTO FINISHED 
END IF 

GOTO checkagain 


KeyBreak: 

COLOR 14 
PRINT "BREAK" 

PRINT 

PRINT "Interrupted by User" 
GOTO StopMotor 

OutOf Bounds : 

' check again 
GetPosition 


IF (ty < HLoLim) 
IF (ty > HHiLim) 
IF (Tz < VLoLim) 
IF (Tz > VHiLim) 
IF (tX < LLoLim) 
IF (tX > LHiLim) 
GOTO checkagain 


THEN GOTO Out 2 
THEN GOTO Out2 
THEN GOTO 0ut2 
THEN GOTO 0ut2 
THEN GOTO Out 2 
THEN GOTO Out2 


Out 2 : 

COLOR 12 

PRINT "BREAK - RIG AT ty; Tz ; tX 
PRINT 

PRINT "Survey Rig Out of Bounds!" 

StopMotor: 

OPEN "O", 3, portl$ 

PRINT #3, "STOP;" 

CLOSE 3 

PRINT "PRESS [R] to RESUME (i.e. glitch on the encoders) any other key to stop" 
a$ = " “ 

WHILE a$ = n “ 
a$ = INKEY$ 

WEND 

IF a$ = "R" OR a$ = "r" THEN GetPosition: GOTO checkagain ELSE GOTO FINISHED 


FINISHED: 


END SUB 

SUB Talk8400 (CMD$) 

' init psi 

IF psi% = 0 THEN CALL IBDEV(0, 5, 0, 13, 1, 0, psi%): PRINT psi% 
IF psi% = -1 THEN PRINT "IEEE488 INIT FALIURE" : STOP 
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REM * send-wait-read to the PSI-8400 

IF echo% = 1 THEN COLOR 10: PRINT , CMD$ ; : COLOR 7 

IF icr% THEN CMD$ = CMD$ + CHR$(13) 

CALL ibwait (psi% , &H8 ) 

CALL IBWRT (psi% , CMD$ ) 

x$ = buffer3$ 

DO 

CALL ibrsp (psi% , spr%) 

LOOP WHILE (spr% AND &H40) = 0 
CALL ibwait (psi%, &H4) 

IF echo% = 1 THEN PRINT 

CALL IBRD (psi% , x$) ' IEEE Read routine 

IF (ibsta% AND EERR) THEN PRINT "IBRD ERROR": STOP 

buffer$ = LEFTS (x$, IBCNTL& ) 

CALL DecodeBuffer (buffers ) 


END SUB 

SUB TalkCompumotor (Hstep&, Vstep&, Lstep&) 

■Set Up and run Utility Program in the Compumotor 3000 


PRINT USING 

" STEP 

: +###### 

+###### +######"; Hstep& ; Vstep&; Lstep& 

I 

OPEN "O", 3 

, portl$ 



; 

PRINT #3, 

"STOP; " 




PRINT #3, 

11 LOAD ; " 



- ..J 

PRINT #3, 

"7000 VELC M ; 

hveloc; “ 11 ; 

LVELOC ; ” “ ; WELOC ; " ; " 


PRINT #3, 

"7001 ACEL M ; 

HAccel; " 

LAccel; " VAccel; ” ; “ 

- 

PRINT #3, 

"7002 MOVE " ; 

Hstepfc; " " ; 

Lstep&; " Vstep&; 


PRINT #3, 

"7003 DONE;" 





PRINT #3, 

PRINT #3, "RUN FROM 7000;" 
CLOSE 3 

END SUB 
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